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Q , Abstract 



This paper is concerned with construction and structural analysis of both cyclic and quasi-cyclic 
codes, particularly LDPC codes. It consists of three parts. The first part shows that a cyclic code given 
by a parity-check matrix in circulant form can be decomposed into descendant cyclic and quasi-cyclic 
codes of various lengths and rates. Some fundamental structural properties of these descendant codes are 
developed, including the characterizations of the roots of the generator polynomial of a cyclic descendant 
code. The second part of the paper shows that cyclic and quasi-cyclic descendant LDPC codes can be 
derived from cyclic finite geometry LDPC codes using the results developed in first part of the paper. 
, This enlarges the repertoire of cyclic LDPC codes. The third part of the paper analyzes the trapping 

i sets of regular LDPC codes whose parity-check matrices satisfy a certain constraint on their rows and 

columns. Several classes of finite geometry and finite field cyclic and quasi-cyclic LDPC codes with 
large minimum weights are shown to have no harmful trapping sets with size smaller than their minimum 
weights. Consequently, their performance error-floors are dominated by their minimum weights. 

L Introduction 

The rapid dominance of LDPC codes [1] in applications requiring error control coding is due to their 
capacity-approaching performance which can be achieved with practically implementable iterative decod- 
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ing algorithms. LDPC codes were first discovered by Gallager in 1962 [1] and then rediscovered in late 
1990's [2], [3]. Ever since their rediscovery, a great deal of research effort has been expended in design, 
construction, structural analysis, efficient encoding and decoding, performance analysis, generalizations 
and applications of LDPC codes. Numerous papers have been published on these subjects. Many LDPC 
codes have been chosen as the standard codes for various next generations of communication systems 
and their applications to digital data storage systems are now being seriously considered and investigated. 

Let GF(q) be a field with q elements. A regular q-ary LDPC code [1] is given by the null space over 
GF(q) of a sparse parity-check matrix H that has constant column weight 7 and constant row weight p. 
Such an LDPC code is said to be (7,p)-regular. If the columns and/or rows of H have varying weights, 
then the null space of H gives an irregular q-ary LDPC code. If H is an array of sparse circulants of 
the same size, then the null space over GF(q) of H gives a g-ary quasi-cyclic (QC)-LDPC code. If H 
consists of a single sparse circulant or a column of sparse circulants of the same size, then the null space 
of H over GF(q) gives a cyclic LDPC code. If q = 2, an LDPC code is said to be binary. 

In almost all of the proposed constructions of LDPC codes, the following constraint on the rows and 
columns of the parity-check matrix H is imposed: no two rows (or two columns) can have more than 
one place where they both have non-zero components. This constraint on the rows and columns of H 
is referred to as the row-column (RC)-constraint. This RC-constraint ensures that the Tanner graph [4] 
of the LDPC code given by the null space of H is free of cycles of length 4 and hence has a girth of 
at least 6 and that the minimum distance of the code is at least 7inin + 1, where 7niin is the minimum 
column weight of H [5], [6]. The distance bound 7inin + 1 is poor for small 7inin and irregular LDPC 
codes, but it is tight for regular LDPC codes whose parity-check matrices have large column weights, 
such as finite geometry LDPC codes [5]-[9], and finite field QC-LDPC codes constructed in [10]-[13] 
and this paper. A parity-check matrix H that satisfies the RC-constraint is called an RC-constrained 
parity-check matrix and the code given by its null space is called an RC-constrained LDPC codes. An 
RC-constrained LDPC code is one-step majority-logic decodable [5], [6]. Furthermore, the RC-constraint 
on the parity-check matrices of LDPC codes allows us to analyze the trapping-set structure [14], [15] 
of RC-constrained LDPC codes which affects their error-floor performances. Analysis of trapping-set 
structure of RC-constrained LDPC codes is a part of investigation in this paper. 

LDPC codes can be classified into two general categories: 1) random or pseudo-random codes that 
are constructed using computer-based algorithms or methods; and 2) algebraic codes that are constructed 
using algebraic or combinatorial tools such as finite fields, finite geometries and experimental designs. 
Codes in these two categories can be classified into two types, codes whose parity-check matrices possess 
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little structure and codes whose parity-check matrices have structures. A code whose parity-check matrix 
possesses no structure beyond being a linear code is problematic in that both encoding and decoding 
implementations become quite complex. A code whose parity-check matrix has structures beyond being a 
linear code is in general more easily implemented. Two desirable structures for hardware implementation 
of encoding and decoding of LDPC codes are cyclic and quasi-cyclic structures. A cyclic LDPC code can 
be efficiently and systematically encoded with a single feedback shift-register with complexity linearly 
proportional to the number of parity-check symbols (or information symbols) [6]. Encoding of a QC- 
LDPC code can also be efficiently implemented but requires multiple shift-registers [16], [17]. It is 
in general more complex than encoding of a cyclic code but still enjoys hnear complexity. However, 
QC-LDPC codes enjoy some advantages in hardware implementation of decoding in terms of wire 
routing [18]. Furthermore, the QC structure allows partially parallel decoding [19] which offers a trade- 
off between decoding complexity and decoding speed, while cyclic structure allows either full parallel 
or serial decoding. In this paper, we show that a cyclic LDPC code can be put in QC form through 
column and row permutations. As a result, a cyclic LDPC code enjoys both encoding and decoding 
implementation advantages. Encoding is carried out in cychc form while decoding is carried out in QC 
form. 

QC-LDPC codes are more commonly studied than cyclic LDPC codes. There are at least a dozen of or 
more methods for constructing QC-LDPC codes, including both algebraic and computer-based methods; 
however, there is only one known class of cychc LDPC codes which are constructed based on finite 
geometries [5]. 

This paper is concerned with constructions and structural analysis of both cyclic and QC codes, 
particularly LDPC codes. It consists of three parts. In the first part, it is shown that a cyclic code given by 
a parity-check matrix in circulant form can be decomposed, through column and row permutations, into 
various cyclic and QC codes, called descendant codes. Some fundamental structures of the descendant 
codes are developed, including the characterization of the roots of the generator polynomial of a cychc 
descendant code. In the second part of the paper, it is shown that RC-constrained cyclic and QC-LDPC 
codes can be derived from the class of cychc finite geometry (FG) LDPC codes based on circulant 
decomposition presented in the first part. Several new families of RC-constrained cychc and QC-LDPC 
codes are presented. The third part of the paper is concerned with trapping sets of RC-constrained regular 
LDPC codes. It is shown that for an RC-constrained (7,/9)-regular LDPC code, its Tanner graph has no 
trapping sets of size smaller than or equal to 7 with numbers of odd-degree check-nodes less than or 
equal to 7. Several classes of cychc and QC-LDPC codes are shown to have large minimum distances (or 
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minimum weights) and no elementary trapping sets [20] with sizes and numbers of degree- 1 check-nodes 
smaller than their minimum weights. 

II. CiRCULANT Decomposition, Cyclic and Quasi-cyclic Codes 

A circulant is a square matrix over a certain field such that every row is the cyclic-shift one place to 
the right (or one place to the left) of the row above it and the first row is the cyclic-shift one place to 
the right (or one place to the left) of the last row. In coding theory, a cychc-shift commonly refers to 
the cycUc-shift one place to the right. Hereafter, by a cycUc-shift, we mean a cycUc-shift one place to 
the right unless expUcitly mentioned otherwise. In this case, every column of a circulant is a downward 
cyclic-shift the column on its left and the first column is the downward cycHc-shift of the last column. 
It is clear that a circulant is uniquely specified (or characterized) by its first row which is called the 
generator of the circulant. The columns and rows of a circulant have the same weight. 



A. Circulant Decomposition 

Let W be an n X n circulant over the field GF(g) where g is a power of a prime. We label the rows 
and columns of W from to n — 1. Let w = {wq, wi, Wn-i) be the generator of W. We denote W 
by ^'(w) = '^{wo,wi,...,Wn-i). Then 



W = ^'(w) 



wq Wi 

Wn-l Wq 
Wi W2 



Wn-2 
Wo 



(1) 



Let '^'(^^(w^) denote the circulant obtained by simultaneously cyclically shifting all the rows of ^'(w) one 
place to the right. Let w^^) denote the n-tuple obtained by cyclic-shifting all the components of w one 
place to the right. Then, it is clear that ^(^^(w) = ^'(w*^^)). Note that 'I'(w) and 'I'(w'^^)) have identical 
set of rows and identical set of columns except that all the columns are cyclically shifted one place to 
the right and all the rows are cychcally shifted upward one place. Therefore, ^'(w) and ^'(w^^)) are 
isomorphic up to cyclic-shift. 

Suppose n can be factored as a product of two positive integers, c and I, such that c 7^ 1 and I ^ 1, 
i.e., n = c-l and c and I are proper factors of n. Let I = {0, 1, 2, • • • , c • Z — 1} be the set of indices (or 
labels) for the rows and columns of the n x n circulant ^(w) given by (1). Define the following index 
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sequences: 



7r(°) = [0,c,2c,--- ,{l-l)c], 
;r=[;rW,7r(°) + l,--- ,7rW+c-l]. 



(2) 
(3) 



Then, tt gives a permutation of the indices in I. Suppose we first permute the columns and then the rows 
of W based on tt. These column and row permutations based on tt result in the following c x c array of 
circulants of size I x I over GF(g): 



$(w) 



*(^)(Wc-l) *(wo) 
*W(W2) *W(W3) 

*W(wi) *W(W2) 



*(Wc_2) 1'(w,,_i) 

*(Wc-3) *(Wc_2) 

*(wo) *(wi) 
*W(Wc-l) *(wo) 



(4) 



where, for < i < c. 



Wi 



*(Wi) = 



Wc+i 



W2c+i 



Wi 



(5) 



(6) 



Each / X / circulant ^(wj) (or 'I'(w,p')) in <I>(w) is called a descendant circulant of ^'(w). Since ^(wi) 
and ^'(w|^'') are isomorphic for < i < c, there are at most c distinct (or non-isomorphic) descendant 
circulants of *(w) in $(w), namely ^'(wo), ^'(wi), • • • ,^'(wc-i). The Z-tuple Wj is called the i-ih 
cyclic section of w. 

Since $(w) is obtained by applying the permutation tt to the columns and rows of the ciruclant ^'(w), 
we write $(w) = 7r(*(w)). Let ir"^ be the inverse permutation of tt. Then ^(w) = 7r~^($(w)). From 
the structure of $(w) displayed by (4), we see that each row of I x I circulants is a right cyclic-shift 
of the row above it, however, when the last circulant on the right is shifted around to the left, all its 
rows are cyclically shifted one place to the right within the circulant. This structure is referred to as the 
doubly cyclic structure which is pertinent to the construction of new cyclic codes, especially new cyclic 
LDPC codes, as will be shown in later sections. From the expression of (4), we see that the descendant 
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circulant ^'(wo) of ^'(w) appears in the array $(w) c times on the main diagonal. For 1 < i < c, the 
descendent circulant ^(wj) appears i times and its shift ^'(^^(wj) (or ^'(Wj^^'')) appears c — i times in 
$(w). ^'(wj) and its shifts appear on an off-diagonal of $(w) starting from the zth position of the first 
row and moving down to the right on a 45° diagonal. When it reaches to the last (rightmost) column of 
$(w), it moves to the left of the next row of $(w) and continues to move down on a 45° diagonal until 
it reaches the last row of $(w). 

Summarizing the above results, we have the following theorem. 

Theorem 1. Given an n x n circulant W = ^(w) over afield with generator w, if n can be properly 
factored, then there is a permutation vr which puts W into an array of circulants of the same size in the 
form of (4). Conversely, if an array ^{w) of circulants of the same size is given in the form (4), then 
there is a permutation 7r~^ which puts the array $(w) into a circulant W with generator w. 

Theorem 1 gives a basis for decomposing a cyclic code into families of cyclic and QC codes or putting 
a group of cyclic codes into a longer cyclic code. 

B. Cyclic and QC Descendants of a Cyclic Code 

In the following, we show that cycUc and QC codes can be derived from a given cycUc code using 
circulant decomposition. The results developed in this section will be used in Section IV to construct 
new cychc and QC-LDPC codes from cyclic FG-LDPC codes. 

Let Cc be an (n,n — r) cyclic code over GV{q) given by the null space of an n x tt, circulant parity-check 
matrix 'H.circ = ^(w) over GF(g) with rank r where w is the generator of the circulant. (For every cyclic 
code, a circulant parity-check matrix Hcirc can always be constructed by using its parity-check vector 
as the generator w of the circulant [15]. This will be reviewed in the next section.) Suppose n can be 
properly factored as the product of two integers, c and I. Then, as shown in Section 11. A, the circulant 
parity-check matrix YLcirc = *(w) of Cc can be decomposed as a c x c array H^c of circulants of size 
Z X Z in the form given by (4) through column and row permutations: 



^'(wo) 



^'(wi) 
*(wo) 

^'(^nWc-l) 



^'(W2) 

^'(wi) 



^(Wc-l) 
*(Wc-3) 



H,c = $(w) 



(7) 



^r(l)(wi) 



*W(W2) 



^r(l)(w3) 



*(wo) 
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where, for < z < c, Wj and *(wj) are given by (5) and (6). Then, the null space of Hgc = $(w) gives 
an (n,n — r) QC code Cqc over GF(q) which is combinatorially equivalent to Cc- We say that {Cc, Cqc} 
form an equivalent pair. Notation-wise, we express Cqc and Cc as Cgc = 7r(Cc) and Cc = 7r~^(Cqc), 
respectively. 

From the array H^c = ^*(w), we can construct new cyclic codes of three different types. These new 
cycUc codes are called cyclic descendant codes (simply descendants) of the cyclic code Cc- The cyclic 
code Cc itself is called the mother code. 

For < i < c, if *(wi) is a nonzero circulant, then the null space over GF(g) of ^'(wi) gives a 
cyclic descendant of Cc, denoted by cf'^ , of length I. This descendant code is referred to as a type- 
1 cyclic descendant of Cc- Since there are at most c distinct non-isomorphic descendant circulants of 
Hcir-c = *(w) in the array H^c = $(w). There are at most c distinct type-1 cychc descendants of Cc- 

From (7), we see that each column of the array Hgc = <I>(w) consists of the circulants in the first 
row of Hgc- For < i < c, each circulant ^'(wj) or its cyclic shift ^''^^^(wj) appears once and only 
once. Since a circulant *I'(wi) and its cyclic shift ^'(^^(wj) differ only in permutation of their rows and 
hence their null spaces are identical. Consequently, the null spaces of all the columns of Hgc = $(w) 
are the same. In fact, the null space of each column of Hgc = $(w) is identical to the null space of the 
following cZ X Z matrix: 

*(wo) 



H 



col 



^'(Wi) 



*(Wc-l) 

For 1 < < c, let ii, i2, . . . , ifc be distinct integers such that < ii, i2, . . . , ifc < c. Let 



H 



col,k 



(8) 



which is a submatrix of ^coi- The null space of ^coi,k gives a cyclic code of length I, denoted by Q ' , 
which is referred to as a type-2 cyclic descendant of the mother cyclic code Cc- 

For 1 < A: < c, let ii, i2, • • • , ifc be a set of distinct integers such that < ii, i2, • • • , ifc < c. Suppose 
we replace the descendant circulants, ^'(wjj, ^'(wjj, . . . , ^(w^) of 'H.circ = ^(w) and all their cyclic 
shifts in the array H^c = $(w^) (see (7)) by zero matrices of size Ixl (if ii = 0, we replace c copies of 
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the circulant, ^'(wo), by c zero matrices). By doing this, we obtain a c x c array tlqc,mask = ^{'^)mask 
of circulants and zero matrices of size I xl. Since the cycUc shift of a zero matrix is also a zero matrix, 
the array ^{w)mask is still in the form given by (4). Then Udrcmask = ^{^)mask = T^~^{^{^)mask) 
gives a new n x n circulant over GF(q) . Let r-mask be the rank of ilcirc,mask = '^{'^)mask- Then the 
null space of Udrcmask = *(w)masfc givcs an (n,n - r^ask) cycUc code c|^^^^ which is referred to as a 
type-3 cyclic descendant of the mother cyclic code Cc- The replacement of a set of circulants in the array 
Hgc = $(w) by a set of zero matrices is called masking [6], [10], [11]. tldrcmask = *(w)masfe and 
Hqcmask = ^{'^)mask Called masked circulant md masked array of Hare = *(w) andH^c = $(w), 
respectively. It is clear that different masking pattern results in a different cyclic descendant code of Cc- 
In Section HI, we will characterize the roots of the generator polynomials of cyclic descendant codes of 
all three types. 

For any pair (s,t) of integers with 1 < s,t < c, let Hqc{s,t) he a s x t subarray of H^c = #(w). 
Since Hqds, t) is an array of circulants, its null space gives a QC code. This QC code is called a QC 
descendant code of Cc (or Cqc). 

C. Cyclic- and QC-LDPC Codes Derived From a Cyclic LDPC Code 

If the circulant parity-check matrix YLcirc = *(w) of Cc is a sparse circulant over GF(g) and satisfies 
the RC-constraint, then the null space of "ti-circ = *(w) gives an RC-constrained cycUc-LDPC code 
over GF(g) . Since the c x c array H^c = ^*(w) is obtained from Hci^c = ^(w) by column and row 
permutations, it also satisfies the RC-constraint. Hence, the null space of Hgc = $(w) gives an RC- 
constrained QC-LDPC code Cqc which is equivalent to the cycHc LDPC code Cc- Since the entire array 
Hqe = $(w) satisfies the RC-constraint, any subarray of H^c = <^*(w) also satisfies the RC-constraint. 
Consequently, all the cyclic descendant codes derived from the cyclic -LDPC code Cc are cyclic-LDPC 
codes, i.e., the null space of the ith descendant circulant ^(wj) (or 'I''^^)(wj)) of 'H.circ = ^'(w) in 
the array Hgc = $(w) gives a cycUc-LDPC code of length the null space of the parity-check matrix 
Hcoi,fe given by (8) gives a cyclic-LDPC code of length I, and the null space of a c x c masked circulant 
= '^{w)mask of Hdrc = ^(w) gives a cyclic-LDPC code of length n. The Tanner graphs of 
the cycUc descendant LDPC codes of Cc have a girth of length at least 6. 

For any pair {s,t) of integers with 1 < s,t < c, let Hgc(s, t) be a s x t subarray of H^c = $(w). Then 
the null space of Hqds, t) gives a QC-LDPC code whose Tanner graph has a girth of at least 6. 

Among the classes of LDPC codes that have been constructed or designed, the only class of LDPC 
codes that are cyclic is the class of finite geometry (FG) LDPC codes [5] whose parity-check matrices 
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are circulants and satisfy the RC-constraint. Cyclic FG-LDPC codes have large minimum distances (or 
weights) and perform well with iterative decoding based on behef propagation. Cychc-LDPC codes 
constructed based on two-dimensional projective geometries have been proved that their Tanner graphs 
do not have trapping sets of sizes smaller than their minimum weights [20]. As a result, their error-floors 
are mainly determined by their minimum weights. Since they have large minimum weights, their error- 
floors are expected to be very low. In Section VII, we will show that the Tanner graphs of the cychc-LDPC 
codes constructed based on two-dimensional Euclidean geometries also do not have trapping sets with 
sizes smaller than their minimum weights. Unfortunately, cyclic FG-LDPC codes form a small class of 
cyclic-LDPC codes. However, using circulant decomposition presented in this section, we can construct 
large classes of cyclic and QC descendant LDPC codes from cyclic FG-LDPC codes, as will be shown 
in Sections IV, V and VI. These cycUc and QC descendant LDPC codes of cycUc FG-LDPC codes also 
have good trapping set structures. 

Construction of QC-EG-LDPC codes through decomposition of a single circulant constructed based on 
hnes of a two-dimensional Euchdean geometry was proposed earlier by Kamiya and Sasaki [9]. In this 
paper, their focus was mainly on construction of high rate QC-LDPC codes and analysis of the ranks of 
their parity-check matrices. In this paper, we propose constructions of both cyclic- and QC-LDPC codes 
through decomposition of a single or multiple circulants constructed based on two and higher dimensional 
Euclidean and projective geometries. We particularly emphasize on construction of cyclic LDPC codes 
and characterization of the roots of their generator polynomials. 

III. Decomposition of Cyclic Codes and Characterization of Their Cyclic 

Descendants 

In this section, we first show that a circulant parity-check matrix of a given cyclic code can be expressed 
as a linear sum of circulants which correspond to the roots of the generator polynomial of the given code. 
From this linear sum of circulants, we then characterize the roots of the generator polynomials of the 
cyclic descendants of the given cycUc code. 

A. Circulant Parity-Check Matrices of Cyclic Codes 

For any positive integer m, let GF(g"^) be an extension field of GV{q). Let Cc be an {n,k) cycUc code 
over GF(g) where ra is a factor of — 1 and (n, q) = 1. Every codeword v = (i;o, ui, • • • , Vn-i) in Cc 
is represented by a polynomial v(X) = + viX + • • • + over GF(g) with degree n — 1 or 

less. The polynomial v(X) is called a code polynomial. An {n,k) cyclic code Cc over GF(g) is uniquely 
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specified by its generator polynomial g(^) = go+giX-\ \- gn^k^iX'"~^~^ + X'"~'' which is a monic 

polynomial of degree n — k over GF(q) and divides — 1 [6], [21]-[24] where 7^ 0. A polynomial 
of degree n — 1 or less over GF(q) is a code polynomial if and only if it is divisible by g(X). Hence, 
every code polynomial v(X) is a multiple of g{X). 

The generator polynomial g{X) of Cc has n — k roots in GF((j'™). The condition (n, g) = 1 ensures 
that all the roots of X" — 1 are distinct elements of GF(g™) and hence all the roots of g,{X) are distinct 
elements of GF(g"*). In the construction of a cychc code, its generator polynomial is often specified by 
its roots. This is the case for BCH and RS codes [6], [21]-[24]. 



where hj G GF(g) for < j < A;, /i^ = 1 and /to 7^ 0. The polynomial h(X) is called the parity-check 
polynomial of C. Let 



Let 



h(X) 



{X- - l)/g(X) 



= hQ + hiX + ■■■ + hkX^ 



(9) 



h(X) 



ho + hiX + ■■■ + hkX'' 



= X''h{X-^) = hk + hk-iX + ■■■ + hoX\ 



(10) 



which is the reciprocal polynomial of h(X). Comparing the coefficients of h(-'^) and h{X), we have 



ho — hk, hi — hk-i 



hk — ho- 



(11) 



Form the following n-tuple over GF(g): 




(12) 



k+l 



where the first k + l components are the coefficients of h(X) and last n — k—1 components are zeros. 



I.e., 



hk+i — hk+2 



hn-l = 0. 



(13) 
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Using the n-tuple h of (12) as the generator, we form the following n x n circulant over GF(g): 



^(h) 



ho hi /i2 
hn-i ho hi 



hn-2 hn-1 
hn-3 hn-2 



hi h2 hs ■■■ hn-1 ho 
In terms of the coefficients of h{X), ^'(h) is given as follows: 

hk hk-i hk-2 ■■■ hi ho 
hk hk-i hk-2 ■■■hi ho 



(14) 



*(h) 





ho 













hk hk-i 



hk hk-i 



hk-i hk-2 hk-3 



ho 












hk 



(15) 



The first n — k rows of Hdrc = 'l'(h) are linearly independent which give the conventional parity-check 
matrix H of the (n,k) cyclic code Cc- The other k rows of Hdrc = 'l'(h) are redundant rows (or linearly 
dependent on the the first n — k rows). Since Hdrc is a redundant expansion of H, the null spaces of H 
and Ucirc give the same cyclic code Cc- The n-tuple h = {hk, hk-i, ■ ■ ■ , /lo, 0, 0, • • • ,0) is commonly 
referred to as the parity-check vector. 

Note that every row (or every column) of the circulant parity-check matrix Hdrc = ^(h) of Cc has a 
zero-span of length n — k — 1 (i.e., n — k — 1 consecutive zeros). It is proved in [25] that this zero-span has 
maximum length and is unique. The maximum zero-spans of different rows of Hdrc start from different 
positions (or different columns). It is shown in [25] that using the parity-check matrix in circulnat form, 
an (n,k) cyclic code C can correct bursts of errors up to the code's burst-correction capability or it can 
correct any burst of erasures of length n — k or less using iterative decoding [15], [25]. Decomposition 
of a burst-error correction cyclic codes gives new burst-error correction cychc descendant codes. 

Suppose that n can be properly factored as the product of two positive integers, c and I. Then ^'(h) can 
be decomposed into acxc array of Ixl circulants in the form given by (4) by applying the permutation 
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TT (defined by (3)) to the columns and rows of *(h), 



$(h) 



^(ho) *(hi) 
*(l)(he-l) ^-(110) 



*(hc-3) ^-(110-2) 



where, for < j < c, 



^r(l)(hi) *«(h2) ••• *W(hi) ^'(ho) 

ht = hk-t for < t < k, 
ht = for t> k. 



(16) 



(17) 
(18) 
(19) 



The null space of $(h) gives a QC code Cqc that is combinatorially equivalent to Cc- 

In code construction, the generator polynomial g(-'^) of an (n,k) cychc code Cc over GF(q) is specified 
by its roots [6], [21]-[24]. Let /3o, A, • • • ,Pn-k-i be the roots of g(X) . Then 



0<j<n-fc 



(20) 



Since g{X)\X"' — 1, — 1) and {n,q) = 1, /3o,/3i,-'' , l^n-k-i are distinct nonzero elements of 

GF(g'"). Let q be a primitive nth root of unity. Then, for < i < n — fc, /3i is a power of a. Since 
a" = 1, (A)" = 1 for < i < n — /u. A polynomial c{X) of degree n — 1 or less over GF(q) is a code 
polynomial if and only if c{X) has Po, Pi, - ■ ■ , Pn-k-i as roots, i.e., c(/3j) = for < ? < n — A;. 

In terms of the roots of g(^), the parity-check matrix of Cc generated by g(-^) is conventionally given 
by the following {n — k) x n matrix over GF(g"*): 



V 







vo 




Vl 




Vn-fe-1 





f3o 
Pi 



PI 



Po 



n-1 



n-l 



on-1 
Pn-k-1 



(21) 



1 Pn-k-l Pn-k-1 ■ 

The rows are linearly independent over GF(g™). An n-tuple over GF(g), c = (co,ci,--- ,c„_i), is a 
codeword in Cc if and only if c • = 0. This is to say that the null space over GF(g) of V gives the 
cyclic code Cc- The null spaces of the circulant parity-check matrix Hdrc and V give the same code Cc- 
The parity-check matrix of Cc in the form of (21) is commonly used for algebraic decoding, such as the 
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Berlekamp-Massey algorithm for decoding BCH and RS codes [6], [21]- [24]. 

In the following, we develop some structural properties of the circulant parity-check matrix Hdrc = 
'J'(h) of Cc- One such structural property is that Hdrc can be expressed in terms of the circulants formed 
by the rows of V. For < i < n — k, let 

v, = (l,A,A',--- ,Pr')- (22) 

be the ith row of V and *(vi) be the n x n circulant over GF(g"*) with as the generator Since 
vo, vi, • • • , Vn-k-i are hnearly independent, the circulants, ^'(vo), *(vi), • • • , *(v„_fc_i), are also hn- 

early independent (i.e., fora^ G GF(g"^) withO <i< n—k, ao*(vo)+ai*(vi)H |-a„_fe_i*(v„_fc_i) ^ 

unless ao = ai = • • • = a„_fe = 0). 
For < i < n — A;, let 

Vi{X) = l + PiX + ^fX^ + ■■■ + ^f-^X^-i (23) 
be the polynomial representation of ith row of V and 

Wi{X) = ^f-i + Pl'-^X + ■■■ + PiX""-^ + (24) 
be the reciprocal of Vj(X). For < z < n — A;, since 

X" - 1 = (X - A)(/3r~^ + pr'^x +■■■+ + ^n-i)^ 

then we have 

vn 1 

Vi(X) = — — = + + • • • + A^"-' + X^-\ (25) 

X — Pi 

It follows from (9), (20), partial-fraction expansion and (25) that the parity-check polynomial h(X) of 
Cc can be expressed as a linear combination of Vi{X)s as follows: 



n x-Pi 

0<i<n-k 

aiiX"" - 1) 



n—k— I 



= Yl ^^v^(^) (26) 
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where for < i < n — k, 

n (/^^-/^i) • (27) 

j=o,j^i ) 

Since /3o, /3i, • " " > Pn-k-i are distinct nonzero elements of GF(g™), all the coefficients, cto, ai, • • • , an-k-i, 
of the Unear sum of (26) are nonzero. 

Summarizing the above results, we have the following theorem. 

Theorem 2. Let Cc be an (n,k) cyclic code over GF(q) generated by g{X) which has the following 
nonzero elements of GF(q™-), /3q,(3i,--- , f3n-k-i> roots. For < i < n — k, let Vj(X) = f3^~^ + 
+ • • • + /3jX"~^ + X"'~^. Then the parity-check polynomial h(X) of Cc can be expressed as a 
linear sum of \q{X),yi{X), ■ ■ ■ , v„_fe_i(X) as follows: 

n—k—l 

h(X)= <^MX), (28) 

where, for < i < n — k, 

n—k—l 
j=0,j^i 

Replacing X in (28) by X~^, multiplying both sides by X'^~^, using (10) and (23), the expression of 
(28) can be put in the following form: 

n—k—l 

Xn-k-i^x)= J] am{X). (30) 
The vector representation of the polynomial X"'~''~^h{X) is 

y^{n-k-l) = (0, 0, . . . , 0, hk, hk-U ■■■ , ho), (31) 

which is the (n— A;— l)th right cyclic-shift of the vector representation h = (%,%_!,•• • ,ho,0,0,--- ,0) 
of the reciprocal polynomial h(X) of the parity-check polynomial h{X) of Cc- Putting (30) in vector 
form, we have 

n—k—l 

h(n-fe-l) = ^.v,. (32) 

1=0 

If we cyclically shift the components of all the vectors in (32) k + 1 places to the right, then we have 

n—k—l 



1=0 
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where 



= (34) 

is the {k + l)th right cyclic-shift of Vj, for < z < n — k. It follows from (33) and (34) that we have 

n—k—l 



h= 5^ Xi^i (35) 



1=0 

where, for < i < n — A;, 

K = '^iPi'"'^- (36) 
Then, it follows from (28), (35) and (36) that we have Theorem 3. 

Theorem 3. For an (n,k) cyclic code Cc over GF(q) whose generator polynomial has elements Pq, Pi, - ■ ■ , Pn-k-i 
of GF( q"^), as roots, then 

rt— fc— 1 

h(X)= J2 ^i^i(^) (37) 



where, for < i < n — k, 



1=0 



n—k—l 



= /3r'~M n • (38) 

The circulant parity-check matrix Hdrc of Cc given by (15) can be expressed as the following linear sum 
of circulants, ^'(vo), *(vi), • • • ,^(v„_fc_i), 

n—k—l 

Hei.c = ^(h)= (39) 

i=0 

where for < i < n - k, Vi = {1, Pi, pf, ■■■ , P^'^). 

The circulants, *(vo), *(vi), . . . , ^'(v„_fe-i), are called the root circulants of the cyclic code Cc- It 
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follows from (9), (10), (23) and (37) that the coefficients of the parity-check polynomial h{X) are: 

n—k—l 

hj = Yl ^i^i^ f"'' < i < A;, (40) 
hj = 0, for k < j <n. (41) 

B. Characterization of Cyclic Descendants of a Cyclic Code 

In the following, we characterize the roots of the generator polynomial of a cyclic descendant of 
an (n,k) cyclic code Cc over GF(g) whose parity-check matrix is given in terms of roots of the form 
given by (21). Consider the circulant *(vj) with Vj = {l,Pi,/3f, - ■ ■ , as the generator. Decompose 

*(vi) into a c X c array of Z x Z circulants. The descendant circulants in the first row of ^'(vi) are 
*(vi,o), *(vi,i), • • • , ^'(vi,c-i) where for < j < c, 

Vi,o = ,Pt'^"), (42) 

Vij = /3/v.,o. (43) 

If follows from (42) and (43) that we have 

The equality of (44) implies that if ^'(vj o) is known, all the descendant circulants, ^'(vjj)'s and ^'(Vj^ j^)'s 
can be constructed from ^'(vj^o) using (44). 

It follows from Theorem 3 that the circulant generated by hj is given as follows: 

n—k—l n—k—l 

^ihj)= Xi^i^ij)= Y ^^^'*(vi,o), (45) 



1=0 



(1) 



where hj, the jih cychc section of h, is given by (17). The null space of *(hj) gives a cyclic code Cj 
over GF(q) of length I, a type-1 descendant of Cc- 

For < 11,12 < n — k, suppose there exists an integer t with < t < c such that jSi^ = a^'^Pi^. In 
this case, since a'^' = a" = 1, we must have (3'^^ = Pf^ . We say that and are equal in cth power. 
Then, it follows from (42) and (44) that Vi^^o = ^i2,o and ^'(vjj o) = ^(vj2,o)- Let m be the number 
of distinct circulants among 'I'(vo,o)! ^(vi,o), • • • , ^(v„_fc_i,o)- Then, we can partition the n — k — 1 
roots, Po,Pi, - ■ ■ , Pn-k-i, into m equal classes in cth power. For < e < m, let 

Qe = {Pe,0, Pe,l, " " " , Pe,re-l} (46) 
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be the eth class of equal roots in cth power where each fi^j in is one of the roots, Po,Pi, - ■ ■ , Pn-k-i, 
and re is the number of equal roots in fig- It is clear that 1 < rg < c. For < / < rg, let 

K,f = {l,l3lf,P!:f,--- (47) 

Since = (51, = ■■■ = we have v^^o = v* ^ = ■ • • = v* and ^{i^l,) = ^{i^l,) = ■■■ = 

*(Ve,re-i)- For < / < re, *(v*j) is one of the circulants *(vo,o), *(vi,o), • • • , *(vn-ik-i,o) in the 
second sum of (45). For < e < m, let 

L = {Xe,o, Xe,i, • • • ) Ag^^e-i} (48) 

be the set of coefficients, Aj, of the circulants, 'I'(v*q), 'I'(v* i), • • • , 'I'(v*^ _^), in the second sum of 
(45). Grouping the identical circulants in the second sum of (45) together and for each e with < e < m, 
using ^'(v* q) to represent the eth group of identical circulants, we have 

n—k—l 
rn—l 

= E<.*(^e,o)> (49) 

e=0 

where 

r-e-l 

From (49), we see that the circulant ^'(hj) with generator hj is a linear sum of the m circulants, 
*('^o,o)' ^(vi,o)) ■ ■ ■ ) ^(Vto-i,o)' where for < e < m, the circulant *(Vg q) is generated by v* q = 
(l,/3^o,/32c , . . . Then, it follows from (21), (28), (37), (38), (42), (49) and (50) that we have 

the following theorem. 

Theorem 4. The generator polynomial gj^^ {X) of the type-1 cyclic descendant code C^p of the cyclic 
mother code Cc given by the null space of the I x I circulant *(hj) has P^ q, < e < m, as a root if 
and only if A* ^ ^ 0. 

Theorem 4 characterizes the roots of the generator polynomial of a type-1 cychc descendant of a given 
cyclic mother code Cg. 

Example 1. Let a be a primitive element of GF(2^^). Consider the binary primitive (2047,2025) BCH 
code whose generator polynomial ^{X) has a, 0^,0:^,0:^ and their conjugates as roots. The length 2047 
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of the code can be factored as a product of c = 89 and I = 23. The 2047 x 2047 circulant parity-check 
matrix Hcirc of this BCH code can be decomposed into an 89 x 89 array $(h) ofcirculants of size 23 x 23 
by column and row permutations tt defined by (3). The null space of each 23 x 23 descendant circulant 
of <l>(h) gives the (23,12) Golay code with generator polynomial I + X + X'-' + + X'^ + X'^ + X^'^ 
[6], which has f3 = a^^ , P'^ , P'^ , and their conjugates as roots. A A 

(2) 

The next theorem characterizes a type-2 cychc descendant Q ' of Cc given by the null space of the 
parity-check matrix flcoi,k- 

Theorem 5. For 1 < k < c, let 11,12, ■ ■ ■ ,ik be a set of distinct integers such that < ii,i2, ■ ■ ■ ,ik < c. 
For 1 < t < k, let gj^^ {X) be the generator polynomial of if-th type-1 cyclic descendant code C^^^ of 

Cc given by the null space of ifth descendant circulant *(wjj ofUdrc = *(w). Then the generator 

(2) (2) 
polynomial gj, (X) of the type-2 cyclic descendant code of Cc given by the null space of the 

parity-check matrix ^col,k of (8) is the least common multiple of g^^^"* (X), g^^ (X), . . . , g|^^ {X), i.e., 

gf{X) = LCM{gf}{X), Q<t<k}. (51) 
The roots of g^^\x) is the union of the roots of gj^\x),gj^\x), . . . ,g\^\x). 

(3) (3) 

Consider the parity-check matrix H^^^ mask ^ type-3 cyclic descendent C^^^j^. of Cc- Express each 

(3) 

row of Hcirc mask ^ polynomial of degree n — 1 or less with the leftmost component as the constant 
term and the rightmost component as the coefficient of the term of degree n — 1. This polynomial is call a 
row polynomial. Find the greatest common divisor h^asjt(X) of all the row polynomials. Let ^^^ski-^) 
be the reciprocal polynomial of h^asfe(X). Then the generator polynomial of C^^^j^. is given by 

g2L(^) = - ^)Mmlkin (52) 

IV. Decomposition of Cyclic Euclidean Geometry LDPC Codes 

In this section, we give constructions of new cychc and QC-LDPC codes by decomposing the circulant 
parity-check matrices of the cyclic Euchdean geometry (EG) LDPC codes. 

A. Cyclic Descendants of Two-Dimensional EG-LDPC Codes 

Consider a two-dimensional Euclidean geometry EG(2,g) over the field GF(g), where g is a power of a 
prime [6], [22], [26]. This geometry consists of points and q{q + l) lines. A point in EG(2,g) is simply 
a two-tuple a = (ao, ai) over GF(g) and the zero two-tuple (0,0) is called the origin. A line in EG(2,g) is 
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simply a one-dimensional subspace, or its coset, of the vector space of all the two-tuples over GF(g). 
A line contains q points. If a point a is on a Une £ in ^G{2,q), we say the Hne C passes through a. Any 
two points in EG(2,g) are connected by one and only one line. For every point a in EG(2,g), there are 
(g + 1) lines that intersect at (or pass through) the point a. These lines are said to form an intersecting 
bundle of Unes at the point a. For each line in EG(2,g), there are q — 1 lines parallel to it. Two parallel 
Hnes do not have any point in common. The q{q + 1) lines in EG(2,g) can be partitioned into {q + 1) 
groups, each group consists of q parallel lines. A group of q parallel hnes is called a parallel bundle. 

The field GF(g^), as an extension field of the ground field GF(g), is a reahzation of EG(2,g). Let a 
be a primitive element of GF(g^). Then, the powers of a, a~°° = 0, = 1, a, a^, • • • , af''^, give 
all the elements of GF(g^) and they represent the q^ points of EG(2,g). The 0-element represents the 
origin of 'EG{2,q). 

Let EG*(2,g) be the subgeometry obtained from 'EG{l,q) by removing the origin and the g + 1 lines 
passing through the origin. This subgeometry consists of — 1 non-origin points and — 1 lines not 
passing through the origin. Each line in EG*(2,g) has only q — 2 lines parallel to it. Hence, each parallel 
bundle of lines in EG*(2,g) consists of g — 1 parallel lines not passing through the origin. Each intersecting 
bundle of lines at a non-origin point consists of q hnes. Let L = {a-'S a^^, • • • , a^"} be a hne in 
EG*(m,g). For < z < g2 _ 1, let a^C = {q;-'^+% • • • , a^"-^'}. Then, a'C is also a hne in 

EG*(2,g) and q;°>C, aC, ■ ■ ■ , ofl^~'^L give all the (f' — 'V lines in EG*(2,g). This structure of lines is 
called cyclic structure [6], [7]. 

Let £ be a line EG*(2,q'). Based on L, we define the following (g^ — l)-tuple over GF(2), 

whose components correspond to the q^ — 1 non-origin points a°, a, a^, • • • , a^^~^ of EG*(2,(7), where 
Vj = 1 if is a point on jC. and Vj = otherwise. It is clear that the weight of vc is q- This (g^ — l)-tuple 
V£ is called the incidence vector of the hne £ [5], [6]. Due to the cychc structure of the hnes in EG*(2,g) 
(i.e., if £ is a line, a£, is also a line), the incidence vector Vq,£ of the line aC is the cychc-shift (one 
place to the right) of the incidence vector V£ of the line jC. 

Let n = q'^ — 1. Form an n x n matrix Heg over GF(2) with the incidence vectors of the n hnes, 
a^C, ajC, • • • , a^~^jC, of EG*(2,q) as rows. Then, H^;^ is an n x n circulant with both column and row 
weights q. Heg can be obtained by using the incidence vector vc of the line C as the generator and 
cychcally shifting ^^c n — 1 times. Since two hnes in EG*(2,gf) have at most one point in common, their 
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incidence vectors have at most one position where they both have 1 -components. Consequently, Heg 
satisfies the RC-constraint and its null space gives a cychc EG-LDPC code Ceg [5], [6], [15] whose 
Tanner graph is free of cycles of length 4 and hence has a girth of at least 6. The RC-constraint on the 

parity-check matrix H^;^ ensures that the minimum weight (or distance) of Ceg is at least q + 1. To find 
the generator polynomial gEciX) of Ceg, we express each row of Heg as a polynomial over GF(2) 
of degree n — 1 or less with leftmost entry as the constant term and rightmost entry as the coefficient of 
X"--^. Let h{X) be greatest common divisor of the row polynomials of H^;^. The reciprocal liEGi^) 
^eg{X) is the parity-check polynomial. Then, the generator polynomial gEG{X) = (X" — l)/h.EG{X). 

For the special case with g = 2'', the rank of H^g is 3* — 1 [27], [6], [9] and the minimum weight 
of Ceg is exactly 2^* -|- 1 [6], [28]. An integer h with < ^ < 2^*, can be expressed in radix-2* 
form as follows: h = cq + ci2*, where < co,ci < 2*. The sum W2s{h) = co + ci is called the 
2*-weight of h. For any non-negative integer I, let /i^'^ be the remainder resulting from dividing 2'/i by 
2^' - 1. Then < h^^^ < 2"^' - 1. The radix-2'' form and 2^* -weight of /i^) are /i^) = cj^ + cfh" and 
W2={h^''^) = Cq'' + cf \ respectively. Then, is root of the generator polynomial gEG{X) of Ceg if 
and only if [6], [28] 

< max Wsaf/i^')) < 2^ (53) 

0<l<s 

The smallest integer that does not satisfy the condition given by (53) is 2* -|- 1. Hence, gEGi^) has the 
following consecutive powers of a, a^, • • • , as roots. 

Constructions of cyclic LDPC codes based on finite geometries, Euclidean and projective, were first 
presented in [5]. In [5], the authors showed that cyclic finite geometry (FG) codes perform very well over 
the AWGN channel with iterative decoding based on behef propagation (IDBP) using the sum-product 
algorithm (SPA) and the decoding of these codes converges very fast. 

Let c and I be two proper factors of n such that n = c-l. Decompose the n x n circulant parity-check 
matrix Heg into a c x c array 7r(H^c') = 7r('I'(v£)) of circulants of size Ixl in the form of (4) through 
column and row permutation vr defined by (3). Note that every row of -k(J1eg), as a {q'^ — 1) x (g^ — 1) 
matrix, still corresponds to a line in EG*(2,g) not passing through the origin of EG(2,g). Since Heg 
satisfies the RC-constraint, each descendant circulant in ttCHeg) also satisfies the RC-constraint. 

Based on the array Tr{tlEG) of circulants, three types of cychc descendant LDPC-codes of the cyclic 
EG-LDPC code Ceg can be constructed. Note that the first row of Heg is not the parity-check vector. 
For q = 2*, the roots of the generator polynomial gEG{X) of Ceg can be determined from (53). Then, 
it follows from Theorems 4 and 5, the roots of the generator polynomials of a type-1 and type-2 cychc 
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descendant codes can be determined. QC-EG-LDPC codes can also be constructed by taking the null 
spaces of subarrays of ttCHeg)- 

For g = 2*, let 2* — 1 = c • I. Let v = {vf),vi, ...,V2=-2) be the incidence vector of a chosen 
line in EG(2,2'') not passing through the origin as the generator of the (2* ~ 1) x (2* — 1) circulant 
Heg = over GF(2). For < i < c, let Vj = (-Uj, Vc+j, f be a cyclic section of v. 

The ranks of Heg = ^(v) and its type-1 circulant descendant ^'(vj) and type-3 circulant descendant 
tlEG,mask = *(v)masfc (masked circulant of Ueg = ^(v)) as defined in Section n.B can be determined 
easily. Let a be a primitive element of GF(2*). Define the following two (2* — 1) x (2'' — 1) matrices over 
GF(2^): V = [a-'^] and V"^ = [a% <i,j <2^ -1. Both V and V"^ are Vandermonde matrices 
[23], [24] and non-singular. Furthermore, VV~^ = I where I is a (2'' — 1) x (2'' — 1) identity matrix. 
Hence, is the inverse of V and vice versa. Then, the matrix 

2=-2 2=-2 2=-2 

is a (2* — 1) X (2* — 1) diagonal matrix over GF(2*) whose ith diagonal element, < i < 2* — 2, equals 

2''-2 

E a^-^Vj. The vector composed of the diagonal elements of H-^^^ is the Fourier transform [23] of the 

3=0 

incidence vector v = {vq,vi^ ...,V2=-2)- ^eg called the Fourier transform of Ueg- ^eg ^eg 
have the same rank. Since H^g a diagonal matrix, its rank, denoted by rankCH-^Q), is equal to the 
number of nonzero diagonal elements in H^^ which is 3* — 1, same as that of H eg- 

Similarly, the rank of the type-3 circulant descendant Ueg, mask = *(v)masik of H^g; = *(v) is 
equal to the number of nonzero diagonal elements of its Fourier transform {flEG,mask)'^ of ^EG,mask- 

To determine the rank of a type-1 descendant circulant ^'(vj) of Ueg = We define V = 

and = [fJ''-'], < i,j < I where j3 = cc"^. The order of P is /. Then, for < z < c, the Fourier 
transform of ^(vj) is 

(*(vi))-^ =V*(v,)V-i 

i-i i-i i-i 
= diag{ E Vjc+i, E P^Vjc+h • • • , E /3^^"^^^^^jc+i) 

j=0 j=Q j=0 

l-l 

is an i X / diagonal matrix over GF(2*) whose ith diagonal element, < i < /, equals E P^''vjc+i- The vec- 

j=0 

tor composed of the diagonal elements of (^'(vj))-^ is the Fourier transform of Vj = {vi, Vc+i, ■ ■ ■ , V(i_i)c+i) 
Then, for < i < c, (^'(vj))-^ is the Fourier transform of the type-1 descendant circulant ^'(vi) of 
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Heg = ^(v). (^'(vi))-^ and ^'(vi) have the same rank. Hence the rank, ranA;(*(vi)), is equal to the 
number of nonzero diagonal elements in (^'(vj))-^. 

To determine the rank of the parity-check matrix tlcoi,k of a type-2 cyclic descendant code given by 
(8). We first find the Fourier transform of each I x / circulant descendant in Hcoi,k- Divide the rows of 
the Fourier transforms of the k descendant circulants in iicoi,k into / groups. Each group 9j, I < j < I, 
consists of the jth rows of the k descendant circulants in tlcol,k- A group is called a nonzero group if 
not all its k rows are zero rows, otherwise called a zero group. Then the rank of tlcoi,k is equal to the 
number of nonzero groups of rows in the Fourier transforms of the k descendant circulants in Hco;,jfc. 

Example 2. Let the two-dimensional Euclidean geometry EG(2, 2^) over GF(2^) be the code construction 
geometry. The field GF(2^'^) is a realization ofEG(2, 2^). Based on the incidence vectors of the 2^^^ — 1 = 
4095 lines not passing the origin of EG(2, 2^), we can construct a 4095 x 4095 RC-constrained circulant 
Heg 'with both column and row weights 64. Any line not passing through the origin of EG(2,2^) can 
be used to construct the generator (the first row) of Heg- The rank of VLeg is 3^ — 1 = 728. The 
null space of Heg gives a (4095,3367) cyclic EG-LDPC code Ceg 'with minimum distance 65. Its error 
performances decoded with 50 iterations of the sum-product algorithm (SPA) [3], [6], [15] and the 
scaled min-sum (MS) algorithm [29] over the binary AWGN channel are shown in Figure 1. We see 
that the error performance of the code decoded with 50 iterations of SPA is slightly better than that 
of 50 iterations of the scaled MS-algorithm. Furthermore, decoding of the code with the MS algorithm 
converges very fast. The performance curves with 5, 10 and 50 iterations of the scaled MS-algorithm 
almost overlap with each other. Also included in Figure 1 is the error performance of the code decoded 
with the soft-reliability based iterative majority-logic decoding (SRBI-MLGD) devised in [30]. We see 
that, at bit-error rate (BER) of 10~^, the SRBI-MLGD performs only 0.6 dB from the scaled MS with 
50 iterations. The SRBI-MLGD requires only integer and binary logical operations with a computational 
complexity much less than that of the SPA and the MS-algorithm. It offers more effective trade-off between 
error-performance and decoding complexity compared to the other reliability-based iterative decoding, 
such as the weighted bit-flipping (WBF) algorithms [5], [6], [15], [31], [32]. 

AA 

Example 3. Consider the 4095 x 4095 circulant Heg constructed in Example 2. Suppose we factor 
4095 as the product ofc = 3 and I = 1365. By column and row permutations, the 4095 x 4095 circulant 

Heg can be decomposed into a 3 x 3 array 7r(H£G') of descendant circulants of size 1365 x 1365 in 
the form given in (4). Let ^q, and ^2 denote the 3 descendant circulants of Heg in the flrst row of 



23 



ttCHeg)- Then 



7I"(Heg) = 



*0 *1 *2 



r/ie descendant circulants and ^2 both have column and row weights 24. The descendant circulant 
has both column and row weights 16. The rank of^i is 600 (the number of nonzero diagonal elements 
of its Fourier transform '^f). Consider the cyclic LDPC code given by the null space of^i. This 
code is a (1365,765) cyclic EG-LDPC code with rate 0.56 and minimum weight at least 17, the column 
weight of *i plus 1. The code is a type-1 cyclic descendant of the cyclic (4095,3367) EG-LDPC code 
given in Example 2. Its generator polynomial has P = a^,P'^, ...,,0^^ consecutive power of P as roots 
where a is a primitive element of GF(2^^). It follows from the BCH bound [6], that the minimum weight 
is again at least 1 7 which agrees with bound of column weight plus one. By extensive computer search, 
we find that C^j^q has no trapping set with size smaller than 1 7 (see Section VI), however, we do find 
a (17,0) trapping set which gives a codeword of weight 17. Therefore, the minimum weight of C^q is 
exactly 17 and the error-floor of this code is dominated by the minimum weight of the code. The error 
performance of the code over the AWGN channel using BPSK signaling decoded with 50 iterations of 
the SPA (or MSA) is shown in Figure 2(a). At the block error rate (BLER) of 10^^, the code performs 
1.6 dB from the sphere packing bound. 
Suppose we use 



as a parity-check matrix. This matrix is a 4095 x 1365 matrix over GF(2) with constant column weight 

64 but two different row weights, 16 and 24. Its rank is 664 and hence it has a large row redundancy 

(2) 

(3431 redundant rows). The null space of'H.coi,3 gives a (1365,701) cyclic-EG-LDPC code C^^q with rate 
0.5135 and minimum distance at least 65. It is a type-2 cyclic descendant of the (4095,3367) cyclic EG- 
LDPC code given in Example 2. The error performances of this code over the AWGN channel decoded 
with 50 iterations of the SPA and the SRBI-MLGD-algorithm are shown in Figure 2(b). This code is 
one-step majority-logic decodable and it can corrects 32 errors with simple one-step ( OS) majority-logic 
decoding (MLGD) [6]. 

Suppose we replace the circulants, ^'2 and its cyclic-shift ^'^\ in Tr(H.EG) by two 1365 x 1365 zero 



^0 




^2 
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matrices O. We obtain the following 3x3 masked array of circulants of size 1365 x 1365; 





*0 


^-1 


o 




o 










O 





The above array is still in the form of (4) with doubly cyclic structure. It is a 4095 x 4095 matrix 
over GF(2) with both column and row weights 40. Applying the inverse permutation tt"^ to the rows and 
columns ofTr{'H.EG)rnask> we obtain an RC-constrained 4095x4:095 circulant H^iG.mask with both column 
and row weights 40. The rank ofiH^Q mask is J 392. The null space ofYlECmask gives a (40,40)-regular 
(4095,2703) cyclic-EG-LDPC code with minimum distance at least 41. It is a type-3 cyclic descendant 
code of the (4095,3367) cyclic EG-LDPC code given in Example 2. The error performances of this code 
over the AWGN channel decoded with 3, 5 and 50 iterations of the SPA is shown in Figure 2(c). A A 

We can factor 4095 as the product of 15 and 273. Setting c = 15 and / = 273, we can decompose 
the 4095 x 4095 circulant H^;^ given in Example 2 into a 15 x 15 array ttCHeg) of circulants of size 
273 X 273. From this array of circulants, we can construct many type- 1,-2 and -3 cyclic descendant LDPC 
codes of the (4095-3367) cychc EG-LDPC code Ceg given by the null space of Heg- 

In this section, we have shown that given a two-dimensional Euchdean geometry, many cychc EG- 
LDPC codes with large minimum weights can be constructed. 

B. Quasi-Cyclic Descendants of Two-Dimensional Cyclic EG-LDPC Codes 

In the previous subsection, we have considered constructions of cyclic descendant LPDC codes of 
cyclic EG-LDPC codes based on two-dimensional Euclidean geometries. In this subsection, we consider 
constructions of QC descendant LDPC codes of cyclic EG-LDPC codes based on two-dimensional 
Euclidean geometry. As pointed out earher that construction of QC descendant EG-LDPC codes based on 
two-dimensional Euclidean geometries was also proposed in [9]. However, the approach to construction 
proposed in this section is different, mathematically simpler and more general than that in [9]. The 
approach in conjunction with masking allows us to construct both high and low rate codes. Furthermore, 
a fundamental theorem on decomposition of a circulant parity-check matrix H^g constructed based 
on a two dimensional Euclidean geometry into an array of circulant permutation matrices (CPMs) is 
proved. This theorem will be generalized for constructing QC -EG-LDPC codes based on high-dimensional 
Euclidean geometries. Therefore, the construction of QC descendant EG-LDPC codes is a generalization 
of that proposed in [9]. 
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In the following, we will present two types of QC descendant EG-LDPC codes. First, we consider the 
RC-constrained c x c array ttCFLeg) of circulants over GF(2) of size of Z x Z constructed in the previous 
subsection where cl = n = — 1 and I > q — 1. For a pair of positive integers, (s,t) with 1 < s,t < c, 
let 7r(H£;G)(s, t) be a s X t subarray of ttCHeg)- This subarray also satisfies the RC-constraint and its 
null space gives a QC descendant LDPC code of the cyclic EG-LDPC code Ceg given by the null 

space of the ra X n circulant Heg- The QC-LDPC code is referred to as a type-1 QC descendant 

code of Ceg- Note that q does not divide — 1. For I > q — 1 and be a factor of — 1, the smallest I 
is q + 1. 

Notice that the transpose of the parity-check matrix tlcoi,k of a type-2 cychc descendant EG-LDPC 
code gives the parity-check matrix Tr{'H.EG){^,k) of a type-1 QC descendant EG-LDPC code. Both 
parity-check matrices tlcoi,k Tr{'H.EG){^,k) have the same rank which is equal to the number of 
nonzero groups of rows in the Fourier transforms of the k circulants in iicoi,k (or the number of nonzero 
groups of columns in the Fourier transforms of the k circulants in TT(HEG){^,k)). 

Example 4. Consider the 3 x 3 array ttCHeg) of circulants of size 1365 x 1365 given in Example 3 
constructed based on the two-dimensional Euclidean geometry EG(2,2^). Set s = 1 and t = S. Take the 
first row [^0 ^2] of ttCHeg) as a 1 x 3 subarray Tr(H. eg) {'^,3) of 7r{ilEG)> i-^-> '^(H.EG)i^,3) = 
[^0 ^2] which is the transpose of the parity-check matrix Heo/,3 of the type-2 cyclic LDPC code given 
in Example 3. 7r(H£;G')(l, 3) is a 1365 x 4095 matrix over GF(2) with constant row weight 64 but two 
different column weights 16 and 20. The null space of this subarray gives a (4095,3431) QC-EG-LDPC 
code, a QC descendant of (4095,3367) cyclic EG-LDPC code given in Example 2. The bit and block 
error performances with 3, 5, and 50 iterations of the SPA are shown in Figure 3. A A 

For a type-1 QC descendant of a cyclic EG-LDPC codes Ceg given by the null space of a {q^ — 1) x 
(g^ — 1) circulant H^;^ constructed based on the 2-dimensional Euclidean geometry EG(2,g), the size of 
each circulant in its parity-check matrix is at least q-\-l. 

Next, we consider type-2 QC descendants of Ceg • Suppose q — l can be factored as a product of two 
integers, b and / with 1 < b,l < q, i.e., q — 1 = bl. Then n = — 1 can be factored as the following 
product: n = {q -\- l){q — 1) = {q -\- l)bl. Let c= {q-\- 1)6. Then, the circulant parity-check matrix Heg 
of the cyclic EG-LDPC code Ceg of length n = q'^ — 1 can be decomposed into an RC-constrained 
(q + 1)6 X (g + 1)6 array 7r(H.EG)cpm of circulants over GF(2) of size I x I. Since 7T{'H.EG)cpm is 
obtained from H^;^ by column and row permutations, the rank of TT(H.EG)cpm is the same as the rank 
of Heg. The following theorem gives a fundamental structure of the array tt (Heg) cpm which allows us 
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to construct a large class of QC-LDPC codes which are QC descendants of the cycUc EG-LDPC code 
Ceg- We will show that each circulant in iT{'H.EG)cpm is either a circulant permutation matrix (CPM) 
or a zero matrix of size I x I. (A CPM is a permutation matrix for which each row is the cyclic-shift of 
the row above it and the first row is the cyclic-shift of the last row.) We call the array TT{ilEG)cpm the 
CPM-decomposition of Heg, where the subscript "CPM" stands for "CPM-decomposition". 

Theorem 6. Let Heg be the {q^ — 1) x (^^ — 1) circulant over GF(2) constructed based on the — I 
lines of the two-dimensional Euclidean geometry EG(2,q) over GF(q) not passing through the origin. 
Suppose q—1 can be factored as a product of two integers, b and I with 1 < b,l < q, i.e., q—l = bl. Let 
c = (q + 1)6. Then, Heg can be decomposed as a (g + 1)6 x (g + 1)6 array TT{ilEG)cpm of circulants 
of size I X Each circulant is either an I x I CPM or an I x I zero matrix (ZM). Each row (or column) 
block of 'K(H.EG)cpm consists of exactly q CPMs and {q + 1)6 — q ZMs. 

Proof: It follows from the definition of the incidence vector of a line in EG*(2,g) that the q^ — 1 
columns of Heg correspond to the g^ — 1 non-origin points, = 1, a, ct^, q'^^^, of EG*(2,g). 
Permuting the columns and rows based on the permutation vr defined by (2) and (3), we decompose the 
circulant Heg into a c x c array Tr{'tlEG)cpm of circulants of size Z x Z in the form of (4). For < j < c, 
consider the jth circulant in the first row block of the array TT{'H.EG)cpm- It follows from the column 
permutation tt that the columns of correspond to the non-origin points, , a'^'^^ , a'^'^'^^ , a^'~^^'^+-'. 
Suppose that *j is neither an Z x Z CPM nor an Z x Z ZM. Then, the first row of must have at least 
two 1-components. Let yi = q;'^"^+-' and y2 = a'^"^"*"-' with < Zi < Z2 < Z, be the points that correspond 
to two positions where the first row of *j have 1-components. Then, 



where A = a^'^"'^)'^ which is a nonzero element in GF(g). Since < Z2 — Zi < Z, A 7^ 1. Let y = r/x -|- z 
be the hne in EG*(2,g) that contains the points (or coimects) yi and y2 where x and z are two hnearly 
independent points in EG*(2,g) and rj G GF(g). Then, 



y2 = Ayi, 



(54) 



yi = ?7ix + z, 



(55) 



y2 = + z- 



It follows from (54) and (55) that we have 



y2 = A771X + Az. 



(56) 
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where Xrji is a nonzero element in GF{q). Since A 7^ 1, the point Az is different from the point z. Equality 
(56) implies that y2 is also a point on the line y' = rjx + Az that is parallel to the hne y = r/x + z. 
However, a point cannot be on two parallel lines. Consequently, the first row of cannot have more 
than one 1 -component and "^j is either a CPM or a zero matrix. 

As a {q^ — 1) X — 1) matrix over GF(2), the first row of Heg (the incidence vector of a line in 
EG*(2,g)) has q one-components. Since 7r(H^;G'),p„i is obtained from H^^ through column and row 
permutations, the first row of 7r{'H.EG)cpm< as a (g^ — 1) x (o^ — 1) matrix over GF(2), also has q one- 
components. Based on the result proved above, these q one-components must distribute in q CPMs in 
the first row block of the array Tr{'H.EG)cpm^ one in each. Consequently, the first row block of the array 
Tr{'H.EG)cpm consists of q CPMs and c — q = {q + l)b — q ZMs of size I x I. Since Tr{llEG)cpm has 
the cychc structure as displayed in (4), every row block of the array is the cychc-shift of the row block 
above it and the first row block is the cyclic-shift of the last row block. This cyclic structure implies 
that every row (or column) block of TT(HEG)cpm has q CPMs and {q + l)b — q ZMs. This proves the 
theorem. ■ 

The array 7r{'tlEG)cpm of CPMs and ZMs of size I x I can be used as the base to construct QC-LDPC 
codes. For any pair of integers, (7,p) with 1 < 7, p < (^-1-1)6, let '!r{'H.EG){l, p)cpm be a 7Xp subarray of 
Tr{'H.EG)cpm- It is an RC-constrained jl x pi matrix over GF(2). Then, the null space of 7r(HEG)(7, p)cpm 
gives a QC-EG-LDPC code CEG,qc{l-, p) of length pi whose Tanner graph has a girth of at least 6. If 
iT{tlEG){l 1 p)cpm has constant column and row weights, then CEG,qc{l, p) is a regular QC-EG-LDPC 
code. Otherwise, 7r(H£;G')(7, p)cpm has multiple column and/or row weights. In this case, the null space 
of 'rr{iiEG){l, p)cpm gives an irregular QC-EG-LDPC code. 

Here we consider a very special subclass of type-2 QC descendant LDPC codes of the two-dimensional 
cychc EG-LDPC code Ceg- The entire array TT{ilEG)cpm is a (g^ — 1) x {o^ — 1) matrix over GF(2) 
with both column and row weights equal to q. The null space of 7T(H.EG)cpm gives a QC-EG-LDPC 
code CEG,qc{{(l + 1)^5 (q + 1)^) of length n = q'^ — 1 with minimum distance q + 1. If q = 2*, then the 
rank of Tr{'H.EG)cpm is 3* — 1 (the rank of TT{'H.EG)cpm is the same as that of Heg)- In this case, the 
null space of 'K{llEG)cpm gives a QC-EG-LDPC code with the following parameters: 

Length: n = 4* - 1, 
Dimension = 4^ — 3*, 
Minimum distance = 2* -|- L 

For a given two-dimensional Euclidean geometry EG{2,q) over GF(gr), the above construction gives a 
family of structurally compatible QC-EG-LDPC codes. 
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Each factor I of q—1 results in a CPM-decomposition of the circulant Heg with CPMs of size I x I. 
A special case of CPM-decomposition of Heg is Z = g — 1. In this case, the CPM-decomposition of 

Heg is a (g + 1) X (g + 1) array TT{'HEG)cpm of CPMs and ZMs of size (q — 1) x (g — 1). Each row (or 
column) block of 7TiilEG)cpm consists of q CPMs and one single ZM. There are a total of g + 1 ZMs in 
Tr{ilEG)cpm- In constructing the circulant Heg^ we can choose a line C such that, after decomposition, 
the q + 1 ZMs in TT{ilEG)cpm lie on its main diagonal. This special case with I = q—1 was first presented 
in [8] as an array of permutation matrices (PMs) of size (g — 1) x — 1) and was later formulated as 
an array of CPMs of size {q — 1) x {q — 1) in [9]. 

Example 5. Consider the 4095 x 4095 circulant Heg over GF(2) constructed based the two-dimensional 
Euclidean geometry EG(2,2^) given in Example 2. Factor 2^^^ — 1 = 4095 as the product of q + 1 = 
2^ + 1 = 65 and q - 1 = 2^ - 1 = 63. Let c = 65 and I = 63. Decompose the 4095 x 4095 circulant 
Heg into a 65 x 65 array 7T(H.EG)cpm of CPMs and ZMs of size 63 x 63. Suppose Heg is constructed 
by choosing a line L not passing through the origin of EG(2,2^) such that, after decomposition ofilEG> 
the 65 ZMs of 'rr{'H.EG)cpm on its main diagonal. The null space of 'K(H.EG)cpm gives (4095,3367) 
QC-EG-LDPC code which is combinatorially equivalent to the (4095,3367) cyclic EG-LDPC code given 
in Example 2. Suppose we choose a 6 x 65 subarray 7r(H£;G')(6, 65)cpm of of ■K{'H.EG)cpm- The null 
space of this subarray gives a (4095,3771 ) code with rate 0.921. The error performance of this code 
with 50 iterations of the SPA is shown in Figure 4. At the BLER of 10~^, the (4095,3771) code performs 
0.75 dB from the sphere packing bound. A A 

Example 6. Continue Example 5. Suppose we factor g — 1 = 63 a* the product of 9 and 7. Set 6 = 9, 

1 = 7 and c=(g + l)6 = 65 x9 = 585. Decompose the 4095 x 4095 circulant Heg given in Example 

2 into a 585 x 585 array ■K{)llEG)cprn of CPMs and ZMs of size 7x7. Choose 7 = 72 and p = 585. 
Take a 12 x 585 subarray 7r(HEG')(72, 585)cpm from 'K(H-EG)cpm- The subarray 7r(H£;G)(72, 585)cpm 
is a 504 X 4095 matrix over GF(2). The null space of this matrix gives a (4095,3591) QC-EG-LDPC 
code with rate 0.877 whose error performance over the AWGN decoded with 50 iterations of the SPA is 
shown in Figure 5. AA 

Example 7. In this example, we construct a long high-rate code and show how close the code performs to 
the Shannon limit. Let the two-dimensional Euclidean geometry EG(2,257) over the prime field GF(257) 
be the code construction geometry. Based on the incidence vectors of the lines in EG(2,257) not passing 
through the origin of the geometry, we construct a 66048 x 66048 circulant Heg with both column and 
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row weights 257. The null space o/Heg gives a cyclic-EG-LDPC code of length of 66048 with minimum 
distance at least 258. 

Set c = q + 1 = 257 + 1 = 258 and I = q-1 = 257 - 1 = 256. Decompose Ueg into a 258 x 258 
array Tr(H. eg) cpm of CPMs and ZMs of size 256 x 256. In this CPM-decomposition, every row and 
every column consists of 257 CPMs and a single ZM. Suppose YIeg is constructed by choosing a line 
not passing through the origin of EG(2,2^) such that the 258 ZMs lie on the main diagonal of the array 

'^(H-EG)cpm- 

Let 7 = 4 and p = 128. Take a 4 x 128 subarray 7r(H£;(3)(4, 128)cpmfrom Tr{'H.EG)cpm> avoiding the 
ZMs on the main diagonal ofTr{'H.EG)cpm- This subarray 7r(HEG)(4, 128)cj,m is a 1024 x 32768 matrix 
with column and row weights 4 and 128, respectively. The null space of tt (Reg) {^^ ^"^8) cpm gives a 
(4,128)-regular (32768,31747) QC-EG-LDPC code with rate 0.969. The error performance of this code 
over the AWGN channel decoded with 50 iterations of the SPA is shown in Figure 6. At the BER of 10^, 
the code performs 0.6 dB from the Shannon limit. A A 

If we select a set of CPMs and their cyclic-shifts in Tr{'H.EG)cpm and replace them by zero matrices 
of size Z X Z, we obtain an array T^(H-EG,mask)cpm of CPMs and ZMs which has the form of (4) with 
doubly cyclic structure. Applying inverse permutation 7r~^ to the rows and columns of 7r{HEG,mask)cpm, 
we obtain a (g^ — 1) x (g^ — 1) masked circulant ilEG,mask over GF(2). The null space ilEG,mask gives 
a cyclic-EG-LDPC code of length — 1. 



C. Masking 

For a pair of two positive integers, (7, p) with 1 < 7, p < g -h 1, let 



7r(H£G)(7,p) 



cpm 
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>7-l,0 B^_i^i ■ ■ • lD^_i^p_i 

be a 7 X p subarray of TT{'H.EG)cpm- A set of CPMs in Tr{'H.EG){'y, p)cpm can be replaced by a set of ZMs. 
This replacement is referred to as masking [6], [8], [10], [11], [15]. Masking results in a sparser matrix 
whose associated Tanner graph has fewer edges and hence fewer short cycles and probably a larger girth 
than that of the associated Tanner graph of the original 7 x p subarray 7r{HEG){l, p)cpm- To carry out 
masking, we first design a low density 7 x p matrix Z(7, p) = [z-ij] over GF(2). Then take the following 
matrix product: Tr{MEG){l, p)cpm = Hi, p) ^ T^(H-EG)il, p)cpm = [zij'^ij], where ZijBij = Bij for 



B. 



(57) 
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Zij = 1 and Zij'Bij = 0{a{q — 1) x (g — 1) zero matrix) for Zij = 0. We call Z{'j,p) the masking 
matrix, Tr(H.EG){'y, p)cpm the base array and tt{'M.eg){i, p)cpm the masked array. Since the base array 
7r(H£;G) (7, p)cpm satisfies the RC-constraint, the masked array tt{M.eg){i, p)cpm. also satisfies the RC- 
constraint, regardless of the masking matrix. Hence, the associated Tanner graph of the masked matrix 
7r(M^G')(7, /5),.p,„ has a girth at least 6. The null space of the masked array 7t{'M.eg){i, p)cpm gives a 
new QC-EG-LDPC code. If both the masking matrix and the base array are regular, the masked array is 
also regular and its null space gives a regular QC-LDPC code. However, if the masking matrix is irregular 
and base array is regular, the masked array is irregular and its null space gives an irregular code. A well 
designed masking matrix results in a good LDPC code. Design and construction of masking matrices for 
constructing binary LDPC codes are discussed in [6], [8], [10], [11]. 

Example 8. In this example, we construct a long irregular QC-EG-LDPC code using the masking 
technique presented above. Consider the 258x258 array'K{'H.EG)cpm ofCPMs andZMs of size 256x256 
constructed in Example 7. Take a 128 x 256 subarray 7r(H£;(3)(128, 256)cpm /''O'w '^{)^EG)cjmi- We use 
this subarray as a base array for masking to construct an irregular code of rate 1/2. Next we construct 
a 128 X 256 masking matrix Z(128, 286) (by computer search) with column and row weight distributions 
close to the following variable-node and check-node degree distributions (node perspective) of a Tanner 
graph optimally designed for an irregular code of rate 1/2 and infinite length ( using density evolution 
[33]): 

X{X) = 0.4410X + 0.3603X2 ^ o.OOlTlX^ + 0.03543X6 + 0.09331X^ + 0.0204^^ 
+0.0048X9 + 0.000353X27 + 0.04292X^9, 

and 

p(X) = 0.00842X7 _^ 0.99023X8 + 0.00135X9. 

where the coefficient of X* represents the percentage of nodes with degree i + 1. The column and row 
weight distributions of the constructed masking matrix Z(128, 256) are given below: 

v{X) = 106X + 105X2 + 35X8 ^ 10X^9, 

c(X) = 10X^ + 118X8, 

where the coefficient X* gives the number of columns (or rows) of Z(128, 256) with weight z + 1. 
Masking the 128 x 256 subarray 7r(HBG)(128, 256)cpm with Z(128,256), we obtain a 128 x 256 
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masked array 7r(MEG)(128, 256)cpm = Z(128,256) (g) 7r(HBG)(128, 256)cpm o/256 x 256 CPUs and 
ZMs. It is a 32768 x 65536 matrix over GF(2) with average column and row weights 3.875 and 7.75, 
respectively. The null space o/ 7r(M£;G)(128, 256)cpm gives an irregular (65536,32768) QC-EG-LDPC 
code. The error performance of this code with 50 iterations of the SPA is shown in Figure 7. We see 
that at a BER of 10~^, the code performs 0.6 dB from the Shannon limit without visible error floor. Also 
include in Figure 7 is the performance of the DVB S-2 standard (64800,32400) LDPC code [34] with a 
BCH outer code. The DVB S-2 LDPC code is an IRA (irregular repeat-accumulated) code [15], [35]. 
The BCH code is a (32400,32208) shortened BCH code with error-correction capability 12. The BCH 
outer code is used to push down the error-floor of the DVB S-2 code. We see that the (65536,32768) 
QC-EG-LDPC code outperforms DVB S-2 code with the BCH outer code. AA 

V. Construction of QC-LDPC Codes Based on Decomposition of Multiple Circulants 
Constructed from High-Dimensional Euclidean Geometries 

In the last subsection, we considered decomposition of the single RC-constrained circulant constructed 
based on the lines of a two-dimensional Euclidean geometry ¥.G{2,q) over a finite field GF(q') not passing 
through the origin of the geometry into a {q + l)bx {q + l)b array of CPMs and ZMs of size I x I where 
b and I are factors of g — 1 and bl = q—\. From this array of CPMs and ZMs, we can construct a family 
of RC-constrained QC-EG-LDPC codes of various lengths and rates and a family of cyclic LDPC codes. 

In this section, we consider decomposition of multiple circulants constructed based on hnes of an 
m-dimensional EucUdean geometry EG(m,g) over the Galois field GF(g) into arrays of CPMs and ZMs 
of size I X I. From these arrays, we can construct a very large array of CPMs and ZMs which forms a 
base array to construct a large family of RC-constrained QC-EG-LDPC codes. 

Consider the m-dimensional Euclidean geometry EG(m,g) over G¥{q). This geometry consists of q"^ 
points and J = — l)/{q — 1) lines. Each line consists of q points. The field GF(g'™) as an 

extension field of the ground field GF(q) is a realization of the geometry EG(m,g) [6], [26]. Let a be 
a primitive element of GV{q^). Then, the powers, a~°° = {),a^ = 0, a, ...,q;^'"~^, represent g"* points 
of EG(m,g). Again, the element a~°° = represents the origin of EG(m,g). Let EG*(m,g) be the sub- 
geometry obtained by removing the origin and the Hne passing through the origin from EG(m,g). This 
sub-geometry consists of g"* — 1 non-origin points and Jq = {q^~^ — 1) (g"* — l)/(g — 1) hnes not passing 
through the origin of EG(m,g). 

Let £ = {a-'^, a-'^, with < ji,j2,---,jq < g"* — 1 be a line in EG*(m,g) consisting of 
the points, a^^,a^'^,...,a^''. For < i < g™ — 1, a^C = {a-^^"*"*, a-^^^*, a-''"'"*} is also a line in 
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EG*{m,q) [6], [7], [15]. The lines C,aC,a'^C, ...,a'i""-^C are (g"* - 1) different Unes in EG*(m,g). 
Since a^"*"^ = l,a'^"''~^C = L. The g"* — 1 hnes, C,aC,a'^C,...^a'f"~'^C, are said to form a cyc/ic 
c/a55, denoted by Qc- The Jq lines in EG*(m,g) can be partitioned into Kq = {q^~^ — l)/{q — l) cyclic 
classes. 

For any line C in EG*(m,g) not passing through the origin, the incidence vector of £ is a (g™ — l)-tuple 
over GF(2) defined as follows: V£ = {vq,vi, ...,Vq^-2), whose components correspond to the — 1 
non-origin points, a° = 0, a, ...,Q!^"~^, of EG*(m,g) , where vj = 1 if is a point on £,, otherwise 
Vj = 0. The weight of the incidence vector of a Une is q. Due to the cyclic structure of the lines in 
EG*(m,q), the incidence vector v(a'+^>C) of the line a^'^^jO is right cyclic-shift of the incidence vector 
v{a'C) for < z < g"' - 1. 

Denote the Kq cychc classes of hnes in EG*(m,q) with Qco:Qci, ■■■,Qcko-i- ^^'■^^ cyclic class 
Qc, of q"^ — 1 lines with < z < Kq, we form a (g"* — 1) x (9™ — 1) circulant H^;^,* with the incidence 
vectors of the lines £j,a£j,a^£i, as columns such that each column is downward cyclic- 

shift of the column on its left and the first column is the downward cyclic-shift of the last column. This 
(g"^ — 1) X (g"^ — 1) ciculant H^;^ ^ satisfies the RC-constraint and has both column and row weights 
equal to g. Let q = where p is a prime. For s > 3 and m > 3, g is very small compared to g™ — 1. 
Therefore, Heg j is a very sparse circulant. 

Form the following (g"*— l)xKo(g"*— 1) matrix overGF(2) with circulants, tlEG,i,^EG,2, ■ ■ ■ , ^eg,Ko 
as submatrices: 

HfiCgc = [Heg,o Heg,1! • • • , H^cft-o-i]- (58) 

This matrix has column and row weights g and qKo, respectively. Since the columns of llEG,qc correspond 
to the lines of EG*{m,q), tlEG,qc satisfies the RC-constraint. Its null space gives an RC-constrained 
QC-EG-LDPC code 

Cqc,m of length Ko^q^ — 1) with minimum distance at least g -|- 1. The subscript 
"m" stands for the dimension of the Euclidean geometry EG{m, g) used for code construction. 

Suppose g — 1 can be factored as a product of b and / with < 6, / < g, i.e., g — 1 = 6/. Then g™ — 1 
can be factored as follows: 

q^^-l = (qm-l + qm-2 ^ _ _ _^ q ^ _ 
= (g"*-^ + g"'-^ + ... + g + 1)61. 



Let 



C = (g™-l + g— 2 + ... + g + 1)6. 



(59) 
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Theorem 7. For < z < Kq, each {q^ — 1) x {(f^ — 1) circulant tlEG,i constructed based on ith cyclic 
class Qhi of lines of the sub-geometry EG*(m,q) can be decomposed into a cx c array TT{1lEG,i)cpm 
of CPMs and MZs of size I x I by applying the n-permutation to both the columns and rows of YiEG,i- 
Each row (column) block of T^{ii-EG,i)cpm consists of q CPMs and c — q ZMs. 

Proof: The proof of this theorem is similar to the proof of Theorem 5. ■ 
Again, we call 7r{'H.EG,i)cpm the CPM-decomposition of HEG,i- Replacing each circulant llEG,i in 
(58) by its CPM-decomposition 7r(HEG,i), we obtain the following c x cKq array of CPMs and ZMs 
of size I X I over GF(2): 

''^{'H.EG,qc)cpm = W(H.EG,o)cpm''^(H.EG,l)cpm ■■■ T^(H.EG,Ko-l)cpm]- (60) 

The array Tr{'H.EG,Qc)cpm is a sparse array with relatively small number of CPMs compared to the 
number of ZMs. It also satisfied the RC-constraint. Its null space gives a QC-EG-LDPC code which is 
combinatorially equivalent to the QC-EG-LDPC code Cgc,m given by the null space of llEG,qc of (58). 
For 1 < 7 < c and 1 < /9 < cKq , take a 7 x p suarray Tr{'H.EG,qc){l, p)cpm from 7r(H.EG,qc)cpm- This 
subarray is 7Z x pi matrix over GF(2). Its null space gives a QC-EG-LDPC code of length pi which is 
referred to as a QC descendant of the QC-EG-LDPC code Cqc,m given by the null space of ilEG,qc of 
(58). The above construction gives a large family of QC descendant LDPC codes of Cgc,m- 

Again, a special case is 6 = 1 and Z = g — 1. In this case, c = (g™~^ + q"^-^ -|- ... -|- g -|- 1) and 
Tr{llEG,qc)cpm is a c X cKq array of CPMs and ZMs of size {q — 1) x {q — 1) over GF(2). 

Consider the c x c subarray Tr{'H.EG,i)cpm of CPMs and ZMs. As stated in Theorem 7, each column 
(or row block) consists of q CPMs and c — q ZMs. Suppose q can be factored as a product e and /, 
i.e, q = ef. We can spht each column block of Tr{'H.EG,i)cpm into e column blocks of the same length 
with the q CPMs evenly distributed into the new e column blocks, each with / CPMs. This column 
splitting operation is referred to column block splitting. In distributing the CPMs into e new column 
blocks, their relative positions are not changed. This column block splitting results in a c x ce array 
NLcoi,i{(i) of CPMs and ZMs of size I x I, each column block consisting of / CPMs and each row block 
consisting of q CPMs. Next, we spht each row block of Mco;,i(e) into e new row blocks of the same 
length with the q CPMs evenly distributed among the e new row blocks, each with / CPMs. This row 
sphtting operation is referred to as the row block splitting. This row block splitting of M.f.oi,i{e) results 
in a ce X ce array ^coi,row,i{^j ^) of CPMs and ZMs of size I x I. The array Mcoi,r-o«),i(e, e) is called 
the e X e expansion of Tr{llEG,i)cpm- Each column block and each row block of M.coi,row,i{G, e) consists 
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of / CPMs. If we replace each c x c subarray '!T{'H.EG,i)cpm in ''^(^EG)cpm given by (60) with its e x e 
expansion M.coi,row,i{e, e), we obtain the following ce x ccKq array: 

MEG,qc = [Mcol ,row,0 

(e, e) Mcoi,row,i{e, e) • • • Mcoi ,rcnju,Ko—l 

(e,e)]. (61) 
Note that M.EG,qc has a much smaller density of CPMs than that of the array Tr{ilEG,qc)cpm- 

Example 9. Let q = 2^. Consider the 3 -dimensional Euclidean geometry EG(3,2^) over GF(2^). This 
geometry has g^ — 1 = 2^^^ — 1 = 511 non-origin points and 4599 lines not passing through the origin 
of the geometry. The 4599 lines not passing through the origin can be partitioned into 9 cyclic classes, 
each consisting of 511 lines. Using the incidence vectors of the lines in these 9 cyclic classes, we can 
form 9 circulants, Heg^Oj Hbg^i, Hbq^s, of sized 511 x 511. Factor 511 as the product of b = 73 
and I = q — 1 = 7. It follows from Theorem 6, each 511 x 511 circulant flEG,i can be decomposed into 
a 73 X 73 array Tr{'H.EG,i)cpm of CPMs and ZMs of size 7x7. Each column (row) block consists of 8 
CPMs and 65 ZMs. Form the following 73 x 657 array of CPMs and ZMs of size 7 x 7: 

T^(H.EG,qc)cpm = W(H.EG,o)cpm '^(H.EG,l)cpm ' ' ' '^(H.EG,8)cpm\- 

This array is a 511 x 4599 matrix with column and row weights 8 and 72, respectively. The null space 
of this matrix gives a (8,72)-regular (4599,4227) QC-EG-LDPC code with rate 0.9191. 

Suppose we factor q = 8 as the product of e = 2 and / = 4. Using column and row block splittings, 
each 73 x 73 array Tr{'H.EG,i)cpm can be expanded into a 146 x 146 array yicol,row,i{'2, 2) of CPMs and 
ZMs of size 7x7, each row and column block consisting of 4 CPMs and 142 ZMs. Suppose we take first 
8 of these 146 x 146 arrays and form the following 146 x 1168 array of CPMs and ZMs of size 7 x 7: 

Meg(8) = [Meo,,^o^,0(2,2) Meoi,,.o^,i(2,2) ••• Mcoi,^<,^,7(2,2)]. 

It is a 1022 x 8176 matrix over GF(2) with column and row weight 4 and 32, respectively. The null 
space of this matrix gives a (4,32)-regular (8176,7 1 56) QC-EG-LDPC code with rate 0.8752. This code is 
actually equivalent to the (4,32)-regular QC-EG-LDPC code adopted by NASA as the standard code for 
LANDSAT high-speed communications and other missions [15], [36] where the bit error rate requirement 
is 10~^^. The error performance of this code decoded with 50 iterations of the SPA and 15 iterations of 
the MSA are shown in Figure 8. We see that there is no visible error-floor down to the BER of 10~^^. 
The estimated error-floor of this code is below the BER of 10~^^. At the BER of 10~^^, it performs only 
1.6 dB from the Shannon limit. A hardware decoder for the NASA code has been built. A A 
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VI. Decomposition of Projective Geometry LDPC Codes 

RC-constrained cyclic LDPC codes can also be constructed based on the incidence vectors of lines of 
finite projective geometries. For detail construction of this class of codes, the readers are referred to [5], 
[6], [15]. In the following, we consider the decomposition of a subclass of cyclic projective geometry 
(PG)-LDPC codes constructed based on the lines of two-dimensional projective geometries over finite 
fields (often called projective planes). 

Consider the 2-dimensional projective geometry PG(2,g) pver GF(g). This geometry has n = + q + l 
points and n = q'^ + q + l lines [6], [15], [22], [26]. Each line contains of q + 1 points. Two Unes can have 
at most one point in common. Let a be a primitive element of GF(g^). Since q^ — 1 = {q—l){q'^ +q+l), n 
is a factor of q^ — 1. The n points of PG(2,q) can be represented by the n elements of {a°, a, • • • , a"^^} 
[5], [6], [15]. The q+1 points on a line are represented by the q + 1 elements in {a^, a, - ■ ■ , a"^^}. Let 
£ be a line in PG(2,g). The incidence vector of this line C is an n-tuple over GF(2) defined as follows: 
V£ = {vq, vi, - ■ ■ , Vn-i) where vj = 1 if is a point on C, otherwise vj = for < j < n. Since C 
consists of q + 1 points, the weight of V£ is g + 1. It is known that the cyclic-shift of V£ is the incidence 
of another Une in PG(2,g) [6], [15]. The incidence vector V£ and its n — 1 cychc-shifts are all different 
and give the incidence vectors of all the n Unes in PG(2,g). 

Form an n X n circulant Upo over GF(2) with vc and its n — 1 cychc-shifts as rows. The columns and 
rows of UpG correspond to the points and lines of PG(2,gi), respectively. Both column and row weights of 
HpG are equal to g + 1. Since two lines in a projective geometry can have at most one point in common, 
their incidence vectors can have at most one place where they both have 1 -components. Hence, Upo 
satisfies the RC-constraint. Therefore, the null space of Upc gives an RC-constrained cyclic-PG-LDPC 
code CpG of length n = q"^ + q + l and minimum distance at least q + 2, whose Tanner graph has a girth 
of at least 6. 

For the special case q = 2^ the rank of Hpc is 3* -M [5], [6], [15], [27] and the cyclic PG-LDPC 
code CpG has the following parameters: 1) Length n = 2^* -|- 2* -|- 1; 2) Dimension n — 3^ — 1; 3) 
Minimum distance > 2* -|- 2. The roots of the generator g{X) of Cpc can be determined and are given 
in [5], [6], [37]. 

Let c and I be two proper factors of n such that n = c-l. Then, through column and row permutation 

vr defined by (2) and (3), the circulant tlpc can be decomposed into an RC-constrained c x c array 
7r(HpG') of circulants of size of ^ x /. The null space of each nonzero I x / circulant in 7r(HpG') gives 
an RC-constrained cyclic PG-LDPC code of length I. For any pair (7,p) of integers with 1 < p < I, 
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the null space of any 7 x p subarray of 7r(Hp(5) gives a QC-PG-LDPC code of length pi. 

Example 10. Let the two-dimensional projective geometry PG{2,2^ ) over GF(2^) be the code construction 
geometry. This geometry has {2^^ — l)/(2'' — 1) = 4161 points and 4161 lines. Each line consists of 65 
points. Based on the lines of PG(2,2^), we can construct an RC-constrained 4161 x 4161 circulant Hpo 
with both column and row weights equal to 65. The null space of this ciruclant gives a (65,65 )-regular 
(4161,3431) cyclic PG-LDPC code with minimum distance at least 66. The error performances of this 
code over the AWGN channel decoded with 5, 10 and 50 iterations of the SPA are shown in Figure 9(a). 
We see that the decoding of this code converges very fast Since 4161 can be factored as the product of 
3 and 1387. Let c = 3 and I = 1387. Then Hpc can be decomposed into a 3 x 3 array 7r(HpG') of 
circulants of size 1387 x 1387 in the form of (4). Let ^q, ^1 and ^2 be the 3 circulants in the first row 
block of tt^Rpg). The column and row weights of the circulant are both 19. The null space of^i 
gives an RC-constrained (1387,720) cyclic-PG-LDPC code with minimum distance at least 20. Its error 
performance over the AWGN channel decoded with 50 iterations of SPA is shown in Figure 9(b). A A 

Note that n = + g + 1 is not divisible by g — 1. The PG-circulant Hpo cannot be decomposed 
into an array of CPMs of size {q — 1) x (q — 1). Decomposition of circulants constructed based on 
projective geometries of dimensions higher than two can be carried out similar to the decomposition of 
high dimensional Euclidean geometries, except for the CPM-decomposition. 

VII. Trapping Sets of RC-Constrained LDPC Codes 

It has been observed for most LDPC codes, decoded with iterative message-passing decoding algorithms 
such as the SPA or the MSA, that as the SNR continues to increase, the error probability suddenly 
drops at a rate much slower than that in the region of low to moderate SNR (or even stops to drop, 
i.e., the error performance curve flattens out). This phenomenon, known as error-floor, may preclude 
LDPC codes from applications requiring very low error rates. High error-floors most commonly occur 
for unstructured random or pseudo-random LDPC codes constructed using computer based methods or 
algorithms. Structured LDPC codes constructed algebraically, such as finite geometry and finite field 
LDPC codes [5]-[13], in general have much lower error-floors. 

Ever since the phenomenon of the error-floors of LDPC codes with iterative decoding became known 
[38], a great deal of research effort has been expended in finding its causes and methods to resolve or 
mitigate the error-floor problem. For the AWGN channel, the error-floor of an LDPC code is mostly 
caused by an undesirable structure, known as trapping set [14], [15], in the Tanner graph of the code 
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based on which the decoding is carried out. 



A. Concepts and Definitions 

Let C be an LDPC code of length n given by the null space of a sparse m x n parity-check matrix 
^ ~ [^ij]' 0<i<m, 0<j<n over GF(2) with m rows and n columns. The Tanner graph [4] Q of 
C is a bipartite graph with two sets of nodes, the variable nodes (VNs) and the check nodes (CNs). The 
VNs, denoted by vq,v\, represent the n code bits of a codeword v = {vq,v\, in the 

code and the CNs, denoted by co,ci, ...,c^_i, represent the m (parity) check-sum constraints that the 
code bits must satisfy (they must be all equal to zero). For convenience, we do not distinguish a "code 
bit" and a "VN", or a "check-sum" and a "CN". We will use the notation Vj for both the jth code bit 
and its corresponding VN and the notation Cj for both the zth check-sum and its corresponding CN. A 
VN Vj is connected to a CN q by an edge if and only if the code bit Vj is contained in the check-sum 
Cj. Basically, the VNs correspond to the n columns of the parity-check matrix H and CNs correspond 
to the m rows of H. The jth VN Vj is connected to the ith CN Cj by an edge if and only if ft, = 1. 
The degree dy. of the VN Vj is defined as the number of CNs connected to Vj and the degree dc^ of 
the CN Ci is defined as the number of VNs connected to the CN The degree d^. of the VN Vj is 
simply equal to the number of 1-entries in the jth column of the parity-check matrix H = [/ij j] and the 
degree d^ of the CN Cj is simply equal to the number of 1-entries in the ith row of H = [hi^j\. For a 
(7,/9)-regular LDPC code, all the VNs have the same degree 7 and all the CNs have the same degree p. 
For an irregular code, its Tanner graph has varying VN degrees and/or varying CN degrees. It is clear 
that the number of edges in the Tanner graph of an LDPC code is equal to the total number of 1-entries 
in the code's parity-check matrix H. 

Figure 10(a) shows that the Tanner graph of a (3,3)-regular (7,3) LDPC code given by the null space 
of the following RC-constrained parity-check matrix: 

10 110 
10 110 

10 110 
H= 10 11 

1 1 1 
1 1 1 
1 1 1 
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Definition 1. Let Q be the Tanner graph of a binary LDPC code C given by the null space ofanmxn 
matrix H over GF(2). For 1 < k < n and < t < m, a (k,t) trapping set [14] is a set T{k,t) of 
K VNs in Q which induce a subgraph of Q with exactly r odd-degree CNs (and an arbitrary number 

of even-degree CNs). An elementary {k,t) trapping set [20] is a trapping set for which all CNs in the 
induced subgraph of the Tanner graph have degree one or degree two, and there are exactly r degree-one 
CNs. 

In an elementary trapping set, every CN of degree 1 is connected to a single VN and every CN of 
degree 2 (if exists) is connected to two VNs. Figures 10(b) and 10(c) shows two subgraphs of the Tanner 
graph of a (3,3)-regular LDPC code shown in Figure 10(a) which are induced by a (3,3) trapping set and 
a (4,4) trapping set, respectively. The (3,3) trapping set consists of 3 VNs, vi, V4 and vq. The subgraph 
induced by this trapping set has 3 CNs of degree 1 and 3 CNs of degree 2. Therefore, this trapping set 
is an elementary trapping set. The (4,4) trapping set consists of 4 VNs, vi, V2, V4 and vq. The subgraph 
induced by this trapping set has 3 CNs of degree 1, one CN of degree 3 and 3 CNs of degree 2. 

Suppose, in transmission of a codeword, an error pattern e with k errors at the locations of the k VNs 
of a {k, t) trapping set occurs. This error pattern will cause r parity-check failures (i.e., the check-sums 
are not equal to zeros, because each of these r check-sums contain an odd number of errors in e). In this 
case, for iterative decoding, another decoding iteration must be carried out to correct the failed check- 
sum. Iterative decoding, such as the SPA and MSA, is very susceptible to trapping sets of a code because 
it works locally in a distributed-processing manner Each CN has a local processor unit to process the 
messages received from the VNs connected to it and each VN has a local processor unit to process the 
messages received from the CVs connected to it. Hopefully, these local processor units through iterations 
and message exchanges collect enough information to make a global optimum decision of the transmitted 
code bits. 

In each decoding iteration, we call a CN a satisfied CN if it satisfies its corresponding check-sum 
constraint (i.e., its corresponding check-sum is equal to zero), otherwise, call it an unsatisfied CN. During 
the decoding process, the decoder undergoes state transitions from one state to another until all the CNs 
satisfy their corresponding check-sum constraints or a predetermined maximum number of iterations is 
reached. The ith state of an iterative decoder is represented by the hard-decision sequence obtained at 
the end of ith iteration. In the process of a decoding iteration, the messages from the satisfied CNs try to 
reinforce the current decoder state, while the messages from the unsatisfied CNs try to change some of 
the bit decisions to satisfy their check-sum constraints. If errors affect the n code bits (or the k, VNs) of a 
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{k, t) trapping set T{k, t), the r odd-degree CNs, each connected to an odd number of VNs in T(k, r), 
will not be satisfied while all other CNs will be satisfied. The decoder will succeed in correcting the 
errors in T(k, r) if the messages coming from the unsatisfied CNs connected to the VNs in T(/{, r) are 
strong enough to overcome the (false or inaccurate) messages coming from the satisfied CNs. However, 
this may not be the case if r is small. As a result, the decoder may not converge to a valid codeword even 
if more decoding iterations are performed and this non-convergence of decoding results in an error-floor. 
In this case, we say that the decoder is trapped. 

For the binary-input AWGN channel, error patterns with small number of errors (or low-weight error- 
patterns) are more probable to occur than error patterns with larger number of errors. Consequently, in 
message-passing decoding algorithms, the most harmful (k, r) trapping sets are usually those with small 
values of k, and r, especially when the value of r is very small compared to that of k. Extensive study 
and simulation results [13], [38]-[68] show that the trapping sets that result in high decoding failure rates 
and contribute significantly to high error-floors are those with small values k and small ratios t/k. We 
call these trapping sets small trapping sets. The trapping sets with large values r relative to values k in 
general result in relatively small decoding failure rates and contribute little to error-floor. From extensive 
computer simulations reported in the literature [14], [38]-[68], it has been observed that most trapping 
sets that exert a strong influence on the error-floor are of the elementary trapping sets and trapping sets 
with t/k< 1. 

Besides small trapping sets and their distributions, undetected errors caused by small minimum weight 
of a code also contribute considerably to the error-floor of the code. If there are no trapping sets with 
size K smaller than the minimum weight of an LDPC code, then the error-floor of the code decoded with 

iterative decoding is dominated by the minimum weight of the code. For r = 0, T(k, 0) is a special 
trapping set with no odd-degree CN. Such a trapping set is induced by an error pattern which is identical 
to a codeword of weight k . When such a trapping set occurs, the decoder converges into an incorrect 
codeword and commits an undetected error. In this case, we say that the decoder is trapped into a fixed 
point. 

The notion of a small trapping set given above is loosely defined. A more quantitative definition of 
small trapping set was given by Ladner and Milenkovic [20]. 

Definition 2. A {k, t) trapping set in the Tanner graph of an LDPC code of length n is said to be small 
if K < ^Jn and r < 4k (i.e. the ratio t/k < 4j. 

Since Richardson introduced the notion of trapping sets and their effect on error-floor in 2003 [14], 
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a great deal of research effort has been expended in analyzing the general structure of trapping sets 
of LDPC codes, developing methods (or algorithms) for finding trapping sets (especially the harmful 
ones), techniques to remove small trapping sets, and devising decoding strategies to remove or reducing 
the degrading effect of harmful trapping sets, [13], [38]-[68]. The research effort expended so far still 
leaves the trapping set problem largely unsolved. However, study and extensive computer simulations 
have shown that among the trapping sets contribute significantly to the error-floor, the harmful ones are 
mostly the small trapping sets, especially the small elementary trapping sets with t/k < 1. 

Constructing (or designing) codes to avoid harmful trapping sets to mitigate error-floor problem is 
a hard combinatorial problem, just hke finding the number of the minimum weight codewords (or the 
weight distribution) of a hnear code. Consequently, to lower the error-floor of an LDPC code caused by 
(small) trapping sets, an alternate approach is taken. A most commonly taken approach is a decoder-based 
strategy to remove or reduce the effect of harmful trapping sets on error-floor. Several such decoder- 
based strategies have been recently proposed [53], [54], [56], [58], [61], [67], [68]. Among them, the 
most effective decoding strategy is the backtracking iterative decoding algorithm recently presented in 
[68]. 

B. An Analysis of Trapping Sets of the RC-Constrained LDPC Codes 

In this section, we present an analysis of trapping set structure of an RC-constrained regular LDPC 
code. The analysis is based on the RC-constraint on the rows and columns of the parity-check matrix 
H and its column weight 7. For such an RC-constrained LDPC code, its minimum weight is at least 
7 + 1. We will show that there is no (k, r) trapping set with k < 7 and r < 7. More precisely, any 
trapping set (k, r) with k VNs, if k < 7, then the number of odd-degree CNs is at least 7 + 1, i.e.,r > 7. 
This is to say that for an RC-constrained (7,/9)-regular LDPC code, there is no harmful trapping set with 
size smaller than 7. Particularly, we will show that an RC-constrained (7,p)-regular LDPC code whose 
parity-check matrix has column weight 7 has no small elementary trapping sets of the type defined by 
Definition 2. Cyclic EG- and PG-LDPC codes given in [5] and their cychc and QC descendants given 
in Sections IV and VI of this paper are RC-constrained LDPC codes and whose parity-check matrices 
have large column weights, hence they don't have harmful small trapping sets with size smaller than 7. 
Besides the FG-LDPC codes, LDPC codes constructed based on finite fields and experimental designs 
in [10]-[13], [69]-[76] are also RC-constrained LDPC codes. 

Let C be a binary (7,p)-regular LDPC code of length n given by the null space of an RC-constrained 
m X n matrix H = [/iij]o<i<m,o<j<n over GF(2) with column and row weights 7 and p, respectively. 
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Let ho, hi, hm-i denote the rows of H, where the ith row hj is given by the following n-tuple over 
GF(2): hi = (/ii,o,/ii,i, •••,/ii,n-i), for < z < m. An n-tuple v = {vo,vi, ...,Vn-i) over GF(2) is a 
codeword in C if and only if v • H-^ = (a zero m-tuple). The condition v • H-^ = gives the following 
m constraints on the bits of the codeword v: 

Cj = V • hj = vohifi + vihi^i + ... + Vn-ihi^n-i = 0, (62) 

for < z < m, where v • hj is the inner product of v and hj. The above m hnear sums of code bits are 
called parity-check-sums (or simply check-sums). The m check-sums of the code bits equal to are the 
constraints that the code bits of any codeword must satisfy. 

For < j < n, if hij = 1, then the jth code bit vj participates (or is contained) in the zth check-sum 
Cj given by (62). In this case, we say that the zth check-sum Cj checks on the jth code bit vj of v (or 
the jth code bit Vj of v is checked by the ith check-sum q). Since H has constant column weight 7 , 
there are 7 check-sums containing (or checking on) the code bit ^jj. Since every row of H has weight p, 
each check-sum Cj checks on p code bits. Since H satisfies the RC-constraint, no two different rows of 
H have more than one position where they both have 1-components. This implies that no two different 
code bits, Vj^ and vj^, are simultaneously checked by two different check-sums, Cj^ and Cj^. 

Suppose a codeword v = {vo,vi, ...,Vn-i) in C is transmitted over the binary-input AWGN channel. 
Let z = {zo,Zi,..., Zn-i) over GF(2) be the hard-decision received vector (or sequence). The jth received 
bit Zj of z is simply an estimate of the jth code bit Vj of the transmitted codeword v. If Zj = Vj for 
< j < n, then z = v; otherwise, z contains transmission errors. Therefore, z is an estimate of the 
transmitted codeword v prior channel decoding. Let 

e = (eo,ei, ...,e„_i), 

= {Z0,Zi,...,Zn-l) -\- {vo,Vi,...,Vn-l), 

= (zo + Vo, Zi + Vi -\ \- Zn-l+Vn-l). 

where, for < j < n, ej = zj + Vj and "+" is modulo-2 addition. If zj / vj, then ej = 1 otherwise 
Cj = 0. Therefore, the positions in e where the components equal to "1" are the erroneous positions. The 
n-tuple e gives the pattern of errors contained in the received sequence z and is called the error pattern 
contained in z [6]. Hence z = v -I- e. 

For any decoding algorithm (soft or hard), the first step is to compute the syndrome of z [6], 

s= (so,si,...,Sm-i) = z-H"'", (63) 
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where, 

Si = z • hj = zohi^o + zihi^i + ... + Zn-ihi^n-i, (64) 

for < z < m, which is called a syndrome-sum of z. If s = 0, then the received bits in z satisfy all the 
m check-sum constraints given by (62) and z is a codeword. In this case, the receiver assumes that z is 
the transmitted codeword and accepts it as the decoded codeword. If s 7^ 0, the received bits in z do not 
satisfy all the m check-sum constraints given by (62) and z is not a codeword. In this case, we say that 

errors in z are being detected and the error pattern is called a detectable error pattern. Then an error 
correction process is initiated. Since z = v + e and v • hj = 0, it follows from (64) that each syndrome- 
sum is actually a Unear sum of a set of error bits contained in the received sequence z, 

Sj = e • hj = eohifi + eihi^i + ... + e„_i/ii,„_i, (65) 

If z is error-free, Sj = q = for < z < m. If z is not error-free but the error pattern e happens to 
be identical to a nonzero codeword in C, all the m syndrome-sums will be equal to 0. In this case, the 
received sequence z contains an undetected error pattern and decoding results in an incorrect decoding. 

Decoding process is initiated (or continues in iterative decoding) only if not all the syndrome-sums are 
equal to zero. 

From (65), we see that a syndrome-sum Sj is equal to "1" if and only if the number of nonzero error 
digits checked by Sj (or participate in the sum Sj) is odd. A syndrome-sum Sj is equal to zero if and only 
if either all the error bits checked by Sj are error-free or the number of nonzero error bits checked by Si 
is even. Let e be an error pattern with k nonzero error bits which cause r nonzero syndrome-sums and 
an arbitrary number of zero syndrome-sums. Construct a subgraph Q{k,t) of the Tanner graph Q of the 
code with a set T(/«, r) of k VNs. These k VNs correspond to the k nonzero error bits in the detectable 
error pattern e and are connected to r CNs which correspond to the r nonzero (failed) syndrome-sums 
and/or some CNs which correspond to zero syndrome-sums but are adjacent to the VNs in T(k, r). In 
this subgraph, the CNs corresponding to the nonzero (failed) syndrome-sums have odd degrees and the 
other CNs have even-degrees. This subgraph Q{k,, t) is said to be induced by the detectable error pattern 
e and the set T(k, r) is a trapping set as defined in Definition 1 . 

A syndrome-sum Si that contains an error bit Cj is said to check on Cj. Since each column of the 
parity-check matrix H has column weight 7, there are 7 syndrome-sums checking on every error bit Sj, 
i.e, every error bit is checked by 7 syndrome-sums (or contained in 7 syndrome-sums). Since each row 
of H has weight p, each syndrome-sum checks on p error bits. Since H satisfies the RC-constraint, no 
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two error bits can be checked simultaneously by two syndrome-sums. 
For < i < m and < j < n, we define the following two index sets: 

K = {j--0<j< n, hij = 1}, (66) 

Mj = {i : < i < m, hij = 1}. (67) 

The indices in Mi are simply the locations of the 1 -components in the zth row hj of H. J\fi is called 
the support of hj. The indices in Mj give the rows of H whose jth components are equal to "1". 
We call Mj the support of jth code bit Vj. Since H satisfies the RC-constraint, it is clear that: 1) for 

< 11,12 < m and i\ / i2, M-i^ and J\fi.^ have at most one index in common; and 2) for < ji, j2 < n 
and ji / 22, Mj^ and Mj^ have at most one index in common. Since H has constant column weight 7 
and constant row weight p, \Mj\ = 7 for < j < n and = p for < i < m. 
For < j < ra, define the following set of rows of H: 

^(j) = : i e Mj}. (68) 

Then it follows from the RC-constraint on the rows of H that .4^-^^ has the following structural properties: 
1) every row h^^ in A^^^ has a 1 -component at the position j; (2) any 1 -component at a position other 
than jth position can appear in at most one row in A^^^; and (3) for < Ji,j2 < n, and ji 7^ j2, A^^^"^ 
and A'^^^^ can have at most one row in common. The rows in A^^"* are said to be orthogonal on the jth 
code bit Vj. 

For < j < n, define the following set of syndrome-sums: 

= {^^ = e • : G ^(^■)}. (69) 

Then, the jth error bit Cj of the error pattern e is checked by (contained in) every syndrome-sum in S^-^^ 
and any error bit other than Cj is checked by at most one syndrome-sum in S^^'^ . Each syndrome-sum in 

S'^-'^ can be expressed as follows: for i G Mj, 

The syndrome-sums in are said to be orthogonal on the error bit Cj and are called orthogonal 
syndrome-sums on Cj. The RC-constraint on the parity-check matrix H (or property-3 of A^^^) ensures 
that any two different orthogonal sets S''-^^^ and S^-''^^ can have at most one syndrome-sum in common. 
Basically, under the RC-constraint, if two rows in H have 1 -components at two different positions, then 
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the two rows must be identical. 

Consider an error pattern e = (eo, ei, en-i) with a single error at jth position with < j < n, i.e., 
ej = 1. For this single error pattern, all the 7 syndrome-sums in S^^^ orthogonal on ej are equal to "1". 
Since ej is only checked by the syndrome-sums in S^^\ all the syndrome-sums in any other orthogonal 
syndrome set are error free and equal to zero. Consequently, the trapping set correspond to this single 
error pattern is a (1,7) trapping set T(l,7) with one VN and 7 CNs of degree 1. Each of these the CNs 
is connected to the VN vj and has degree one. It is clear that T(l,7) is an elementary trapping set. If 
7 > 4, it is not a small trapping set of the type defined by Definition 2. 

Next, we consider an error pattern e with two errors at positions, ji and j2, i.e., ej^ = ej^ = 1. Then 
all the 7 syndrome-sums in the orthogonal set S^^'^^ check on ej^ and all the 7 syndrome-sums in the 
orthogonal set S^^^^ check on ej^. If S^^'^^ and S^^^ are disjoint, then all the 7 syndrome-sums in 
and all the 7 syndrome-sums in are equal to "1". In this case, the trapping set corresponds to the 
error pattern e with double errors is a (2,27) elementary trapping set T(2, 27) with 2 VNs and 27 CNs 
of degree 1. If S^^^^ and S^^'^^ are not disjoint, then they have exactly one common syndrome-sum which 
checks on both ej^ and ej^ and hence this common syndrome-sum is equal to zero. In this case, the 
trapping set corresponds to the error-pattern e with double errors is a (2, 2(7 — 1)) elementary trapping 
set with 2 VNs, 2(7 — 1) CNs of degree 1 and one CN of degree 2. For 7 > 5, it follows from Definition 
2 that for either case, the trapping is not small. For 7 > 2, the number of odd-degree CNs is greater than 
7. The above analysis shows that the trapping set corresponding to an error pattern with two errors has 
at least 2(7 — 1) CNs of odd degrees. 

Consider an error pattern e with three errors at the positions, ji, j2 and J3. The trapping set corresponds 
to this error pattern has several possible configurations depending the locations of the three errors. The 
first possible configuration is such that the three errors are checked separately by three mutually disjoint 
orthogonal sets, S^^'\ S^^^^ and S'-^^l In this case, the trapping set corresponding to the error pattern e 
is a (3,37) elementary trapping set T(3, 37) with 3 VNs and 37 CNs of degree 1, no CN with even- 
degree. The second possible configuration is such that two orthogonal sets, say S^^^ and S^^^\ have a 
common syndrome-sum and the third orthogonal set S^^^^ is mutually disjoint with S^^'^^ and S^^\ In 
this case, error bits, ej^ and ej^, are jointly checked by a common syndrome-sum in and and 
the error bit ej^ is checked only by the syndrome-sums in S^^\ Then, the trapping set corresponding to 
this triple-error pattern e is a (3, 87 — 2) elementary trapping set with three VNs, 87 — 2 CNs of degree 
1 and one CN of degree 2. The third possible configuration is such that all three errors are checked by 
a syndrome-sum which is contained in all three orthogonal sets, S^^^\ S^^^^ and S^^^K In this case, all 
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the 7 syndrome-sums in each of the orthogonal sets, S^^\ S^^^\ and S^^'^\ are nonzero. The common 
syndrome-sum in these three orthogonal sets contains the three errors, ej^ , Cj^ , ej^ , and all the other 
syndrome sums contain only one of these three errors. Consequently, the trapping set corresponding to 
the error pattern e is a (3, 87 — 2) tapping set with 37 — 8 CNs of degree 1 and one CN of degree 3 (no 
CN with even degree). The fourth possible configuration of a trapping set corresponding to a triple error 
pattern e is such that all three errors ej^ , ej^ and ej^ are checked by the syndrome-sum common to two 
orthogonal sets, say S^^^ '> and S^^^\ and one error, say ej^ is checked by S^^^^ alone. In this case, the 
trapping set corresponding to this triple error pattern e is a (3,87 — 1) trapping set with 87 — 2 CNs of 
degree 1, one CN of degree 3 and one CN with degree-2. The fifth possible configuration is that S'^-'^^ and 
have a common syndrome sum checking on ej^ and ej^, S^^^^ and S^^^^ have a common syndrome 
sum checking on Cj^ and ej^, and S^^^^ and S^^^") are disjoint. For this conjuration, the trapping set is 
a (3,87 — 4) elementary trapping set, with 87 — 4 CNs of degree- 1 and two CN of degree-2. The sixth 
possible configuration is such that the pair of errors, {ej^, ej^), is checked by the common syndrome-sum 
m and S^^'\ the pair (ej-^,ej^) is checked by the common syndrome-sum in S(-'i) and 5(^3), and 
the pair (e,, .ej^) is checked by S^^^^ and S^^'^\ Corresponding to this configuration, the trapping set is 
a (8, 8(7 — 2)) elementary trapping set with 8(7 — 2) CNs of degree 1 and 3 CNs of degree 2. Consider 
the 6 possible configurations of three errors, the trapping set with minimum number of odd-degree CNs 
is the sixth configuration. In this case, the number of CNs of odd-degree (degree 1) is at least 8(7 — 2). 
If 7 > 8, the number of odd-degree CNs in a trapping set correspond to a triple-error pattern is greater 
than or at least equal 7. For 7 > 4, the trapping sets corresponding to the first 4 configurations are not 
small trapping sets of the types defined by Definition 2. If 7 > 6, then the trapping sets corresponding 
to all 6 configurations are not small trapping sets of the types defined by Definition 2. Among all the 
6 possible configurations of errors, the sixth one has the smallest number of CNs of odd degree. For 
this configuration, the 3 errors are uniformly distributed in pairs among the three sets of orthogonal 
syndrome-sums. This maximizes the number of CNs of degree-2 and hence minimizes the number of 
degree- 1. 

For 1 < t < 7 and < ji,j2,---,jt < 7> consider an error pattern of t errors at the positions, 
ji,j2, ■■■,jt- For large t and 7, to analyze all the possible configurations of trapping sets with t VNs is 
very hard if not impossible. However, a lower bound on the minimum number of odd-degree CNs can 
be derived. A configuration that results in a minimum number of odd-degree CNs is such for < r, 

s < t, every pair [ej^,ej^) of errors is checked by the common syndrome-sum in S'(>) and S^^'l This 
configuration actually maximizes the number of CNs with even degrees (all degree 2). If this configuration 
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exists, it results in a {t, t{'j —{t — 1))) trapping set T{t, t{'j — {t— 1))) with t(7 — {t— 1)) CNs of degree 

1 and (t — 1)! CNs of degree 2. It is an elementary trapping set. Any other configurations of t errors 
would results in a trapping set with larger number of odd-degree CNs with multiple odd-degrees. For 
7 > t, a trapping set with t VNs in the Tanner graph of an RC-constrained (7,p)-regular LDPC code has 
at least t(7 — {t — 1)) CNs of odd-degrees. For t < 7, the number of odd-degree CNs is greater than 7. 
For t = 7, the number of odd-degree CNs is at least 7. Based on Definition 2, if t < ^/n and 7 > t + 3, 
there is no trapping set with size smaller than 7 — 3. If the ratio t/k requirement for a small trapping 
set is t/k < 1, then there is no trapping set with size smaller than 7. 

The above analysis shows that the structure, the sizes and the distribution of trapping sets of the 
Tanner graphs of RC-constrained LDPC codes very much depend on the column weights and orthogonal 
sets of rows of the parity-check matrices of the codes. Basically, for a (7,p)-regular LDPC code whose 
parity-check matrix has constant column weight 7, the RC-constraint on the rows and columns of the 
parity-check matrix ensures that: 1) the minimum weight of the code is at least 7 + 1; 2) the girth of 
the code's Tanner graph is at least 6; 3) there is no trapping set of size smaller than 7 + 1 with number 
of odd-degree CNs smaller than 7, (i.e., a trapping set with number of VNs less than 7 + 1 must have 
at least 7 CNs of odd-degrees connected to it); and 4) no trapping set of the type defined by Definition 

2 with size smaller than 7 — 3. Due to these structural properties, RC-constrained (7,p)-regular LDPC 
codes with large 7 in general have much lower error-floors than the unstructured LDPC codes constructed 
using computer-based method, and furthermore decoding of these codes with iterative message-passing 
algorithms converges very fast, as demonstrated by extensive simulation results given in [5]-[13], [15], 
[16], [69]-[71]. 

Summarizing the above results, we have the following two theorems for trapping sets of an RC- 
constrained LDPC code. 

Theorem 8. For an RC-constrained (^,p)-regular LDPC code with 7 > 1, its Tanner graph contains no 
trapping set of size smaller than 7 + 1 for which the number of odd-degree CNs is smaller than 7. 

Theorem 9. Let be a positive integer such that 7 > 3. For an RC-constrained p)-regular LDPC 
code, its Tanner graph contains no trapping set of size smaller than 7 — 3 for which the number of 
odd-degree CNs is smaller than 4(7 — 3). 

The above results are derived based on only the RC-constraint on the rows and column of a parity- 
check matrix and its constant column weight 7 but not on its row weight. Hence, the results apply to 
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LDPC codes whose parity-check matrix has constant column weight but may have multiple row weights. 

It is important to note that a trapping set induced by an error pattern does not necessarily prevent 
decoding to converge unless the error pattern is uncorrectable to the decoder Only the trapping sets 
induced by the error patterns that are uncorrectable to the decoder may prevent decoding to converge (or 
fail) and cause an error-floor in the code's error performance. For long codes, it is hard (or very much 
time consuming) to identify the configurations of those trapping sets which do trap the decoder and 
prevent decoding to converge. However, extensive simulations in many pubhshed hteratures did observe 
that in general, trapping sets of small size (relative to minimum weight oJrmn of the code) with small 
number of odd-degree CNs are the harmful ones. When an error pattern induces such a small trapping 
set, the number of failed CNs is so small such that the messages generated by these failed CNs may 
not strong enough to overcome the messages coming from the satisfied CNs to make the changes of the 
erroneous VNs checked by the failed CNs to satisfy their check-sum constraints. As a result, iteration 
continues. However, for a trapping set with large number of odd-degree CNs compared to its number of 
VNs, the messages generated by the failed CNs would be strong enough to overcome the messages from 
the satisfied CNs to make appropriate changes of some code bits checked by all the CNs in such a way 
that all the check-sums are zero and decoding converges. If a trapping set is induced by an undetectable 
error pattern, it is of the form T(k, 0), in which the k erroneous VNs form a codeword of weight n. 
When this happens, the resultant syndrome of a hard-decision decoded sequence is zero. In this case, 
decoding stops and the decoder converges to an incorrect codeword. This results in an undetected error. 
If the minimum weight Wmin of an LDPC code is small, trapping sets induced by uncorrectable error 
patterns that correspond to the minimum weight codewords may contribute significantly to the error- 
floor in the performance of the code. Based on the above observation and reasoning, a code, in general, 
has a lower error-floor if it does not have small trapping sets (relative to the minimum weight) with 
small numbers of odd-degree CNs and its minimum weight oJmin is reasonable large. If a code does 
not have trapping sets with size smaller than its minimum weight Wmim then the error-floor of the code 
is dominated by its minimum weight a;iniin i-C, dominated by the trapping sets, T(a;inin, 0)s, induced 
by the error patterns that are identical to the minimum weight codewords. In the next two sections, we 
will show several classes of RC-constrained LDPC codes do have large minimum weights and do not 
have harmful trapping sets with size smaller than their minimum weights. Hence, their error-floors are 
dominated by their minimum weights. 

Since not all the trapping sets defined in Definitions 1 and 2 trap the decoder and prevent decoding 
to converge, the term ,trapping set, is actually misnamed. 
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VIII. Trapping Sets of Cyclic FG-LDPC Codes and Their Cyclic and QC Descendants 

Since cyclic FG-LDPC codes and their cyclic and QC descendants are RC-constrained LDPC codes, 
they have the trapping set structure presented in the last section. In the following, we consider the trapping 
set structures of two special subclasses of cyclic FG-LDPC codes. These two subclasses of FG codes 
have large minimum distances and no small trapping sets. They can be decoded with various decoding 
algorithms ranging from hard-decision, reliability-based to pure soft-decision iterative decoding, such as 
the OSMLGD, the bit-flipping (BF), the weighted BF, the binary message-passing (SRBI-MLGD), the 
min-sum and the sum-product algorithms, to provide a wide range of effective trade-offs between error 
performance and decoding complexity. The first subclass of cycUc FG-LDPC codes is the class of cyclic 
EG-LDPC codes constructed based on the two-dimensional Euclidean geometries presented in Section 
IV. Here, we consider the cyclic EG-LDPC code Ceg (or its QC equivalent Cecqc) constructed based on 
the two-dimensional Euclidean geometry EG(2,2'^) over the field GF(2''). The parity-check matrix H^;^ 
of the code is a (2^* — 1) x (2^* — 1) circulant over GF(2) (or ((7 + 1) x (g + 1) array of (g — 1) x (g — 1) 
CPMs over GF(2)) whose rows are the incidence vectors of the fines in EG(2,2*) not passing through the 
origin of the geometry. The column and row weights of this circular parity-check matrix YLeg are both 
2*. Its rank is 3* — 1. As shown in Section IV. B the null space of gives a (4* — 1,4* — cyclic 
EG-LDPC code Ceg with minimum distance exactly 2"^ -|- 1. With the OSMLGD, this code is capable 
of correcting 2*~^ or fewer random errors (or 2* or fewer erasures). 

Recall that the columns and rows of Heg » as a (2^^* — 1) x (2^^ — 1) matrix over GF(2), correspond 
to the 2^^ — 1 non-origin points and 2^^ — 1 lines (not passing the origin) of EG(2,2''), respectively. The 
symbols of a codeword v = {vq,vi, . . . ,V22s-2) in Ceg correspond to the 2^*^ — 1 non-origin points 
of EG(2,2'') and therefore correspond to the columns of Hgc . Since any two points in EG(2,2*) are 
coimected by a fine, any two code symbols are checked by a row in Heg . Consequently, for any two 
error symbols, e^^ and e^^, in an error pattern e, the two sets of syndrome-sums, and S^^, orthogonal 
on and e^^ have (exactly) one syndrome-sum in common. 

It follows from the trapping set analysis given in the last section, any trapping set corresponding to 
an error pattern with 2* or fewer random errors will induce a subgraph of the Taimer graph of the code 
which contains at least 2* CNs of odd degrees. This is to say that code has no {k,t) trapping set of 
size K smaller than 2^* + 1 with the number r of odd-degree CNs smaller than 2*. This implies that for 
K < 2^ , there is no (k ,r) trapping set with the ratio t/k < 1. If the (commonly used) requirements of 
small value of k and t/k < 1 are used to define a small trapping set, then the cycfic EG-LDPC code 
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Ceg has no trapping set of size smaller than 2^ + 1 (the minimum weight of the code). For k < 2*~^, 
the number r of odd-degree CNs is greater than 2^~^. Since the code is capable of correcting 2''"^ or 
fewer errors with the OSMLGD, all the trapping sets of size equal to or smaller than 2*~^ are un-harmful 
(i.e., they do not prevent decoding to converge or converge to an incorrect codeword) if the OSMLGD is 
performed before each new decoding iteration. Since the length of the code is n = 4^* — 1, the square root 
of n, ^Jn = \/A? — \ K, 2*. For 1 < k < 2^ — 3, it follows from the tapping set analysis given in the last 
section that for a («;,r) trapping set, the number r of CNs of odd-degree is at least k(2* — (k— 1)) > 4k. 
Then, it follows from Definition 2 that the cyclic EG-LDPC code Ceg has no small trapping set of the 
type defined by Definition 2 with size smaller than 2* — 3. 

Summarizing all the results developed in the last and this sections, we have the following parameters 
for the cyclic EG-LDPC code Ceg constructed based on two-dimensional Euclidean geometry EG(2,2*) 
over GF(2'*) has the structure parameters: 1) length 4* — 1; 2) dimension 4* — 3*; 3) minimum distance 
2^ + 1; and 4) no trapping set of size less than 2* + 1 or (2* — 3) with number of odd-degree CNs less 
than 2* (or less than 4 x 2*). 

In fact, there are many trapping sets of size greater than 2* with number of odd-degree CNs much 
greater than 2*. As an example, we consider an error pattern e with 2* -|- 1 errors at the positions, 
jo,ii, • • • ,i2»-i,i2«- Suppose the errors positions joiii) ■ ■ ■ ,i2«-i> correspond to the 2* points o?°, o?^, 
. . ., o?'^'-^ of a fine L in EG(2,2'') not passing through the origin. The position _72« is any other arbitrary 
position and it corresponds to the point o?'^" . In this case, there is a single syndrome-sum contains 2* 
errors at the positions iojjij ■ ■ ■ )i2«-i> and this syndrome-sum equals zero. Since in a finite geometry, 
any two points are connected by a line. Then any error at the position in the set {j'oj Jij • • • 5^2= i} and 
the error at the position j2= are contained in at most one syndrome-sum and they are the only errors 
in sum. (Note that the rows of the parity-check matrix Heg correspond only to the lines not passing 
through the origin.) Consequently, this syndrome-sum is equal to zero. Recall that each position between 
and n — 1 is checked by 2* rows of Heg. Therefore, for each position ji, < i < 2*, there are at 
least 2^ — 2 syndrome-sums contain only one error in the error pattern e. As a result, the trapping set 
induced by the error pattern e consists of 2'' -|- 1 VNs, at least 2'' (2* — 2) CNs of degree- 1, at most 2* 
CNs of degree-2 and one CN of degree 2*. If s > 3, the number of degree- 1 CNs is much larger than 
the number of VNs in the trapping set. This error pattern is correctable with the OSMLGD. 

Using the geometric structures, configurations of some trapping sets of an EG-LDPC code may 
be analyzed. Consider another case. Let e be an error pattern with 2^ + 2 errors at the positions, 
JO)jir'' ) j2^-i! j2»! j2'>+i- Suppose the errors positions jo,ji,---,j2^-i, correspond to the 2* points 
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, CK-'S . . . , a^'^"-^ of a line C in EG(2,2*) not passing through the origin. The positions and j2<'+i 
are two arbitrary positions which correspond to the points a^'^" and a^^'+^. Assume that a^'^' and 0-^^"+^ 
are not on the same Une. Then each point on C may pair with either point a^'^' or a^^^+i appearing on 
a line. Based on this, we can readily see that the trapping set induced by the error pattern e has at least 
2^(2* - 3) CNs of degree- 1, at most 2''+^ CNs of degree-2 and one CN with degree 2^ If points a^^' 
and a^^'+^ he on the same line, then the number of CNs with degree- 1 is at least 2'*(2* "2). In either 
case, for s > 3, the number of odd-degree CNs is much greater than the number of VNs of the trapping 
set. 

Now we consider a more general case. For < r < 2* — 2, consider an error pattern e with 2^ + r 
errors positions at the positions, jo,ji,---, i2»-i, ^2= , • • • , i2«+r-i- Again, we assume that the 2* positions 
jo,ii, ■ ■ ■ ,j2^-i, correspond to the 2* points a^°,a^'^, . . . ,a^^'-'^ of a line £ in EG(2,2*) not passing 
through the origin. Following the same analysis given above, we can easily show that the trapping set 
induced by this error pattern with 2* + r errors consists of at least 2^(2'' — r) CNs of degree- 1 and at 
most r2^ CN's of degree-2. Since r < 2* — 2, the number of degree- 1 CNs is much larger than the 
number of VNs. For the case r = 0, the trapping set induced by the error pattern e whose error locations 
corresponding to the 2* points of a line not passing through the origin of the geometry has exactly 
2^(2* — 1) CNs of degree- 1 and one CN of degree 2*. Since there are 2^* — 1 hnes not passing through 
the origin (the rows of the parity-check matrix Heg are the incidence vectors of these lines), there are 
2^* — 1 such trapping sets of size 2*. For such a trapping set, the number of CNs of degree- 1 is 2'' — 1 
times larger than the number of VNs. Error patterns corresponding to these trapping sets are correctable 
with the OSMLGD. 

Example 11. Consider the (63,37) cyclic EG-LDPC code constructed based on the two-dimensional 
Euclidean geometry EG(2,2^) over GF(2^). The parity-check matrix of this code is a 63 circulant over 
GF(2) with both column and row weights 8. The minimum weight of this code is 9. The code is capable 
of correcting 4 or fewer errors with OSMLGD. By computer search, we have found all the trapping sets 
induced by error patterns with 3 up to 22 errors. Table 1 gives a partial list of the found trapping sets. 
From the Table 1, we see that for k < 9, the number r of odd-degree CNs associated to every trapping 
set is greater than k . For k = 9, there are (9,0) trapping sets which correspond to minimum weight 
codewords of the code. The square root ~ 8. From Table 1, we see that for k < 6, the number r 
of odd-degree CNs associated with each trapping set of size k smaller than 6 is greater than 4k. Then, 
it follows from Definition 2 that the Tanner graph of the code does not contain small trapping set with 



51 



size K < 6 of the type defined by Definition 2. In decoding of the (63,37) cyclic EG-LDPC code using 50 
iterations of the SPA, none of the trapping sets with size smaller than 9 prevents decoding to converge 
(or trap the decoder) and the error patterns corresponding to these trapping sets are all correctable. The 
trapping sets ^(0, 0), T(10,0), T(11,0), T(12,0) and T(14, 0) result in undetected error (incorrect 
decoding). The error performance of the (63,37) cyclic EG-LDPC code is shown in Figure 11. 

Suppose we consider the (255, 175) cyclic EG-LDPC code constructed based on the two-dimensional 
Euclidean geometry EG(2,2^) over GF(2^). This code has minimum weight 17. Extensive computer search 
found no trapping set of size smaller than 17 which prevents decoding to converge or cause decoding 
failure. We found some large trapping sets with very large numbers of odd-degree CNs but are not harmful. 
These trapping sets are: r(16, 102), r(18, 110), r(21, 102), r(30, 120), r(29, 120), and r(33, 130). 
All but T(30, 120) have r > An. Therefore, only the trapping set T(30, 120) is a small trapping set by 
Definition 2. AA 

Example 12. Consider the (4095,3367) Cyclic-EG-LDPC code with minimum weight 65 constructed 
based on the 2-dimensional EG(2,2^) over GF(2^) given in Example 2. The parity-check matrix of this 
code has column weight 64. The Tanner graph of this code has no trapping set of size smaller than 
64 with number of odd-degree CNs smaller than 64. Note that \/4095 64. It follows from Definition 
2 that the code has no trapping set with size smaller than 61. As shown in figure 1, decoding of this 
code with either the SPA or the SMA converges very fast. Consider the (1365,765) cyclic descendant 
of the (4095,3367) cyclic EG-LDPC code given in Example 3. The parity-check matrix of this code is 
a 1365 X 1365 circulant with both column and row weights 16. For this code, any trapping set of size 
smaller than 17 has at least 16 odd-degree CNs associated with it. Note that \/l365 > 17. Based on 
Definition 2, it has no trapping set with size smaller than 13. AA 

Next, we consider the trapping set structure of a cyclic PG-LDPC code Cpc constructed based on 
the 2-dimensional projective geometry PG(2,g) over GF(q) with q = 2'^. The parity-check matrix of this 
code is a (g^ + g + 1) x (g^ + g + 1) circulant Upo over GF(2) with both column and row weights 
equal to g + 1. The null space of Upo gives an RC-constrained cycHc PG-LDPC code Cpo of length 
n = -\- q -\- 1 and minimum weight at least g + 2, whose Taimer graph has a girth of at least 6. Since 
the H PG satisfies the RC-constraint and its column weight is + 1, it follows from the analysis given 
in Section VII that Cpc has no trapping set T(k, r) for which both the size k and the number r of 
odd-degree CNs smaller than q + l. The square root of the length of the code is ^/n « q. For k < q — 2, 
it easy to check that the number r of odd-degree CNs of a trapping set T(k, r) is greater than An. Then, 
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CpG has no trapping set of the type defined by Definition 2 with size smaller than q — 2. The results on 
trapping sets of the cyclic PG-LDPC code are exactly the same obtained in [20] derived in a different 
approach. Our derivation of the results are simply based on the RC-constraint on the parity-check matrix 
which is much simpler and less mathematical. 

For q = 2*', the cychc PG-LDPC code Cpq has the following structural parameters: 1) length n = 
2^s + 2* + 1; 2) dimension n — 3* — 1; 3) minimum weight at least 2* + 2; 4) girth at least 6; 5) no 
trapping set of size less than 2* -|- 2 with number of odd-degree CNs less than 2^ -|- 1; and 6) no trapping 
set of the type defined by Definition 2 with size less than 2* — 2. 

IX. Other RC-Constrained LDPC Codes and their Trapping Sets 

Besides EG- and PG-LDPC codes, there are other classes of structured RC-constrained LDPC codes. 
These classes of codes are either constructed based on finite fields [11]-[13], [69]-[71] or experimental 
designs [72]-[76]. Codes in most of these classes are QC-LDPC codes. Since the parity-check matrices 
of the codes in these classes satisfy the RC-constraint, their trapping sets have the structure as described 
in Section VII. B. The constructions based on finite fields given in [11]-[13], [69]-[71] are of the same 
nature and they give several large classes of RC-constrained QC-LDPC codes. Among them, several 
subclasses have large minimum weights. In this section, we choose the first class of QC-LDPC codes 
given in [11] for illustration of their trapping set structure. 

Consider the first construction of QC-LDPC codes given in [11]. Let a be a primitive element of the 
Galois field GF(g) Then, a"'^ = 0, q;° = 1, a, • • • , q;«~2 give all the elements of GF(g). Let Crs be the 
cyclic (q' — l,2,g — 2) RS code over GF(g) with two information symbols whose generator polynomial 
g(X) has a, a^, • • • , a^~^ as roots. Then, for < z < g — 1, the two {q — l)-tuples over GF(g), 

and 

\i = (a^a^••• ,a'), 

are two nonzero codewords in Crs with weight q — 1. Note that ui, • • • , Uq_2 are cyclic-shifts of Uq. For 
z = 0, vo = (1, 1, • ■ ■ ) !)• The subscript "rs" of Crs stands for "Reed-Solomon". 

For < z < g — 1, Uj — Vo is a codeword in Crs with weight q — 2 (minimum weight). Form the 
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following {q — l)x{q— 1) matrix over GF(q) with uq — Vq, Ui — Vq, • • • , Uq-2 — Vq as rows: 



Wo 




" aO-1 


a — 1 




1 


Wi 




ai-^ - 1 


aO-1 ■ 


■ ai-^ - 


1 


Wq_2 




a-1 


■ 


■ aO-1 





This matrix W, ,, is the matrix (with rows permuted) given by Eq. (4) in [11] for the construction of the 
first class of QC-LDPC codes. Every row (or column) of W^.^ consists q — 2 distinct nonzero elements 
and one 0-element of GF(g). The q — 1 zero entries of W^s lie on its main diagonal. Therefore, both 
column and row weights of W^s are q — 2. This matrix satisfies the following constraint on the Hamming 
distance between two rows [11]: for < z, j < q—l, j and < c, Z < g — 1, the Hamming distance 
between the two (g — l)-tuples over GF(g), Q;"^Wj and a'wj, is at least q — 2, (i.e., a'^Wj and a'"Wj 
differ in at least q — 2 places). This constraint on the rows of matrix W^s is called the row-distance 
(RD)-constraint and W^s is called an RD-constrained matrix. 

Let P he a {q- 1) X {q- 1) CPM whose top row is given by the {q - l)-tuple (010 • • • 0) over GF(2) 
where the components are labeled from to g — 2 and the single 1 -component is located at the 1st 
position. Then P consists of the {q — l)-tuple (010 • • • 0) and its g — 2 right cyclic shifts as rows. For 
1 < i < g, let P* = P X P X ■ • • X P be the product of P with itself i times, called the ith power of P. 
Then, P' is also a (g — 1) x (g — 1) CPM whose top row has a single 1 -component at the ith position. 
For z = g - 1, P«-^ = the {q-l)x {q- 1) identity matrix. Let P° = P*"^ = I^-i. Then the set 
V = {P°, P, P2, • • • , P*~2} of CPMs forms a cychc group of order q — 1 under matrix multiplication 
over GF(2) with pi-^-^ as the multiplicative inverse of P* and P° as the identity element. 

For < z < — 1, we represent the nonzero element a* of GF(q) by the {q — 1) x {q— 1) CPM P*. 
This matrix representation is referred to as the {q — l)-fold binary matrix dispersion (or simply binary 
matrix dispersion) of a*. Since there are q — I nonzero elements in GF(g) and there are exactly q — I 
different CPMs over GF(2) of size {q — I) x {q — 1), there is a one-to-one correspondence between a 
nonzero element of GF(q) and a CPM of size {q — 1) x (q — 1). Therefore, each nonzero element of 
GF(q) is uniquely represented by a CPM of size {q — 1) x {q — 1). For a nonzero element 6 in GF(q), 
we use the notation B((5) to denote its binary matrix dispersion. If (5 = a*, then B((5) = P*. For the 
0-element of GF(g), its binary matrix dispersion is defined as the {q — 1) x {q — 1) ZM, denote P~°°. 

Dispersing each nonzero entry of W^s into a {q — 1) x {q — 1) CPM over GF(2) and each O-entry 
into a (g — 1) X (g — 1) ZM, we obtain the following {q—l)x{q — 1) array of CPMs and/or ZMs over 
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GF(2) of size {q-l)x{q- 1): 



Bn 



Bg_2 B| 



Bi 





B 



q-2 



Bi Bo 



Bn 



(72) 



where Bj = B{a^ — 1) for < j < g — 1. H,,, is called the binary {q — l)-foId array dispersion of W^-^ 
(or simply binary array dispersion of W^s). This array has {q — 1) ZMs which lie on its main diagonal. 
It is a (g — 1)^ X {q— 1)^ matrix over GF(2) with both column and row weights equal to g — 2. Based 
on the RD-constraint on the rows of W^s and the binary CPM matrix dispersions of the entries of W^s, 
it was proved in [10], [11], [69], [71] that H^^, as a (g - 1)^ x (g - 1)^ matrix over GF(2), satisfies the 
RC-constraint. Hence, its associated Taimer graph has a girth of at least 6. The RD-constrained matrix 
Wrs used for constructing the RC-constrained array Urs of CPMs is called the base matrix for array 
dispersion. 

For any pair (7,p) of integers 7 and p with I p < q, let Hrs(7,p) be a 7 x p subarray of H^s. 
Hrs(7, p) is a 7(g — 1) x p{q — 1) matrix over GF(2) which also satisfies the RC-constraint. The null 
space of H-rsil, p) gives a QC-LDPC code Crs,qc of length p{q — 1) with rate at least {p — ^)/ p, whose 
Taimer graph has a girth of at least 6. If Hrs(7, p) does not contain any of the ZMs of H^s, then H^s has 
constant column weight 7 and constant row weight p. In this case, Crs,gc is a (7, p)-regular QC-LDPC 
code. If Hrs(7, p) contains ZM(s) of H^^, it has two different column weights, 7 — 1 and 7, and/or two 
different row weights, p — 1 and p. In this case, the null space of Ursi'j, p) gives a near-regular binary 
QC-LDPC code. 

For a given finite field GF(g), the above construction gives a family of structurally compatible RC- 
constrained QC-LDPC codes. Consequently, the construction gives a large class of binary QC-LDPC 
codes. Since their parity-check matrices satisfy the RC-constraint, they have the same trapping set structure 
presented in VII. B. 

A very special case is the QC-LDPC code Crs,qcj given by the null space of the full array H^.^ with 
q = 2^. For this case, the length of the code is n = (2* — 1)^ and its minimum weight is at least 2* — 1. 
Using the technique presented in [13], we find that the rank of H^s is 



rank(H.rs) = 3^ — 3. 



(73) 



(The derivation of the expression of (73) is given in a separate paper.) Since the column weight of H,. 
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is 2* — 2, it follows from the analysis of trapping set structure of an RC-constrained LDPC code given 
in Vn that for k < 2* — 2, Crs,qcj has no trapping set T(k, r) of size k, with number of odd-degree 
smaller than 2* — 2. Note that ^/n = 2* — 1. Then, for k < 2* — 5, there is no trapping set T{k, t) of the 
type defined by Definition 2 with number of odd-degree CNs smaller than 4k. That is to say that there 
is no trapping set with size smaller than 2* — 5. 

Summarizing the above results, the QC-LDPC code Crs.qcj given by the full array H^^ of (72) for 
g = 2* has the following parameters: 1) length n = (2^ - 1)^; 2) dimension (2* - 1)^ -3^ + 3; 3) 
minimum weight at least 2* — 1; 4) any trapping set T(k, t) with k < 2* — 2 must have more than 2^ — 2 
CNs of odd-degrees; and 5) no trapping sets of the type defined by definition 2 with size smaller than 
2* - 5. 

Example 13. Let GF(2^) be the field for code construction. Based on this field, we can construct a 
31 X 31 array H^s ofCPMs and ZMs of size 31 x 31. H^s is a 961 x 961 matrix over GF(2) with both 
column and row weights 30. The null space o/H^s gives a (30,30)-regular (961,721) QC-LDPC code 
with minimum distance at least 31. This code is the code given in Example 1 of [11]. For this code, any 
trapping set T{k, t) with k < 30 must have more than 30 CNs of odd-degrees. The code has no trapping 
sets of the type defined by Definition 2 with size smaller than 27. None of the trapping sets with size 
smaller than 31 traps the decoder. The error performances of this code with 5, 10 and 50 iterations are 
shown in Figure 12. AA 

Besides the class of RD-constrained base matrices given above, several other classes of RD-constrained 
base matrices for constructing RC-constrained arrays of CPMs have been proposed in [10]-[13]. Based on 
these arrays of CPMs, several large classes of RC-constrained QC-LDPC codes have been constructed. 
Codes in these classes perform well with iterative decoding using either the SPA or MSA. In the following, 
we describe another method for constructing a large class of RD-constrained base matrices for array 
dispersions to construct RC-constrained QC-LDPC codes. This method is based on a class of Latin 
squares over finite fields and is proposed in [13]. 

An array is called a Latin square of order n if each row and each column contains every element of 
a set of n elements exactly once [77]. Latin squares can be constructed from finite fields. Consider the 
field GF(g). Let a be a primitive element of GF(g) and rj be nay nonzero element of GF{q). Form the 
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following q X q matrix over GF(q'): 



Wls = 



ar] — a 



,0 



a^r] — a 
ar] — a 



a* ^77 — a° ^77 — a 



a °°ri — a° a °°77 — a 



(74) 



Then, V\f is is a Latin square of order q over GV{q) . Every element of GF(g) appears in a row and a 
column once and only once. In [13], it was proved that satisfies the RD-constraint. Binary array 

dispersion of V^ls gives a. q x q array H^^ of CPMs and ZMs of size {q — 1) x {q — 1). Each row or 
column of Uls contains one only one ZM. Hi,5 is a q{q — 1) x q{q — 1) matrix over GF(2) with both 
column and row weights g — 1. 

For any pair ('y,p) of positive integers with 1 < 7, p < let £[^5(7, be a 7 x p subarray of 
Hl5. 11ls{i,p) is a j{q — 1) x p{q — 1) matrix matrix over GF(2). If 11ls{i,p) does not contain any 
ZM of Hi5, then HisiliP)-, as a j{q — 1) x p{q — 1) matrix, has column and row weights 7 and p, 
respectively. The null space of J1ls{i,p) gives an RC-constrained (7,p)-regular QC-LDPC code CLS,qc 
of length p{q — 1). The code has the trapping set structure as described in Section VII. 

For q = 2^, the QC-LDPC code CLs,qc,f given by the null space of the full array Uls has the following 
parameters [13]: 

Length: n = 2^(2* - 1), 

Number of parity symbols: n — A; = 3* — 1, 

Minimum distance dmin > 2* + 2. 
It follows from the trapping set analysis given in Section Vn, any trapping set corresponding to an 
error pattern with 2"^ — 1 or fewer random errors will induce a subgraph of the Tanner graph of the code 
which contains at least 2* — 1 CNs of odd degrees. If the requirements of small value of k and t/k < 1 
are used to define a small trapping set, then the QC-LDPC code CLs,qcj has no trapping set of size 
smaller than 2^ — 1. 

Since the length of the code is n = 2*(2* — 1), the square root of n, ^Jn « 2*'. For 1 < k < 2* — 4, the 
number r of CNs of odd-degrees is at least k(2* — 1 — (k — 1)) > 4k. Then, it follows from Definition 2 
that the QC-LDPC code CLs,qcj has no trapping set of the type defined by Definition 2 with size smaller 
than 2* - 4. 



Example 14. The code constructed based on the Latin square of order 32 over GF(2^) is an RC- 
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constrained (992,750) QC-LDPC code with minimum weight at least 34. Extensive computer search 
found no trapping sets with size smaller than 34 that trap the decoder Two trapping sets T(36, 0) are 
found. This says that the minimum weight of the code is 36. Also found are 1595 T(40, 0) trapping sets. 
Since there are no harmful trapping sets with sizes smaller than the minimum weight, the error-floor of 
the code is dominated by the minimum weight of the code which is 36. The error performances of this 
code over the AWGN channel with 50 iterations of the SPA and the MSA are shown in Figure 13. We 
see that there is no visible error-floor down to the BER of 10~^^. At the BLER of 10~^ (decoded with a 
min-sum FPGA decoder), the code performs 1.1 dB from the sphere packing bound. A A 

X. Conclusion and Remarks 

In this paper, we have shown that cyclic and quasi-cycUc descendant codes can be derived from a 
known cyclic code through decomposition of its parity-check matrix in circulant form using column and 
row permutations. We have analyzed some structural properties of descendant cyclic codes of a cyclic 
code, particularly in characterization of the roots of their generator polynomials. By decomposition of 
cyclic finite geometry LDPC codes, we are able to enlarge the repertoire of cyclic finite geometry LDPC 
codes and to construct new quasi-cychc LDPC codes. The cychc and quasi-cychc structures allow the 
implementation of encoding of LDPC codes with simple shift registers with linear complexity. These 
structures also simplify the hardware implementation of LDPC decoders. Quasi-cychc structure simplifies 
wire routing of an LDPC decoder and allows partial parallel decoding that offers a trade-off between 
decoding complexity and decoding delay. We have shown that a cychc LDPC code can be put in quasi- 
cyclic form through column and row permutations and vice versa. In encoding, we use its cyclic form 
and in decoding, we use its quasi-cyclic form. This allows us to have both advantages in encoding 
and decoding implementations. In this paper, we have also analyzed the trapping set structure of LDPC 
codes whose parity-check matrices satisfy the RC-constraint. We have shown that several classes of finite 
geometry and finite field LDPC codes don't have trapping sets with sizes smaller than the minimum 
weights of the codes. The codes in these classes have large minimum weights. Consequently, codes in 
these classes have very low error-floors which are pertinent to some communication and storage systems 
where very low error-rates are required. 

Finally, we would hke to point out that there are two large classes of structured LDPC codes [8], 
[69] which satisfy the RC-constraint but are not quasi cyclic. The class of LDPC codes given in [8] 
was constructed based on finite geometry decomposition and the class of LDPC codes given in [69] was 
constructed based on Reed-Solomon codes with two information symbols. These two classes of codes 
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have large minimum distances. It follows from our trapping set analysis, they don't have trapping sets 
of sizes smaller than their minimum distances. 
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TABLE I 

A PARTIAL LIST OF TRAPPING SETS OF THE (63,37) CYCLIC EG-LDPC CODES 
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Fig. 1. Bit error performances of the binary (4095,3367) cyclic EG-LDPC code given in Example 1 decoded with the SPA 
and the scaled MSA. 




Fig. 2. (a) The error performances of the binary (1365,765) cyclic EG-LDPC code given in Example 3 decoded with 50 
iterations of the SPA and the MSA. 
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Fig. 2. (b) The error performances of the binary (1365,701) cyclic EG-LDPC code given in Example 3 decoded with the MSA 
and the SRBI-MLGD-algorithm. 




Fig. 2. (c) The bit and block error performances of the binary (4095,2703) cyclic EG-LDPC code given in Example 3. 
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Fig. 3. The bit and block error performances of the binary (4095,3431) QC EG-LDPC code given in Example 4. 
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Fig. 4. The bit and block error performance of the binary (4095,3771) QC-LDPC code given in Example 5. 



66 



10"^ 
10"^ 
10^ 



10"^ 
10"^ 
10"' 
10"' 





□ QC LDPC(4095,3591 ), SPA, BER 
-□- QC LDPC(4095,3591), SPA, BLER 
Sphere Packing Bourid 


1 " I 
i b T 




* 1 

? \ 

1 » 1 




- ' 1 
\ 1 
! q 
' 1 
1 1 




! 1 
1 

' 1 


: 

; ; 




□ 



012345678 



Fig. 5. The bit and block error performance of tlie binary (4095,3591) QC-LDPC code given in Example 6. 




Fig. 6. The bit and block error performances of the binary (32768,31747) QC-LDPC code given in Example 7. 
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Fig. 7. The error performances of the binary (65536,32768) QC-LDPC code and the DVB S-2 standard code given in Example 
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Fig. 8. The error performances of the binary (8176,7156) QC-LDPC code given in Example 9. 
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Fig. 9. (a) The error performances of the (4161,3431) cyclic PG-LDPC code given in Example 10 decoded with various number 
of iterations of the SPA. 
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Fig. 9. (b) The error performances of the binary (1387,720) cyclic LDPC code given in Example 10. 
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Fig. 10. (c) A (4,4) trapping set. 
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Fig. 11. The bit and block error performances of the (63,37) cyclic EG-LDPC code given in Example 11. 
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Fig. 12. The bit error performance of tiie binary (961,721) QC-LDPC code given in Example 13 decoded with 5, 10 and 50 
iterations of the SPA. 




Fig. 13. The bit and block error performances of the binary (992,750) QC-LDPC code given in Example 14. 



